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3.1  INTRODUCTIOSr 


This  chapter  reviews  the  mathematical  tools  and  techniques  required  to 
solve  differential  equations.  Stucfy  of  these  operations  is  a  prerequisite  for 
courses  in  aircraft  flying  qualities  and  linear  control  systems  taught  at  the 
USAF  Test  Pilot  School.  Only  analysis  and  solution  techniques  vdiich  have 
direct  application  for  work  at  the  School  will  be  covered. 

Many  systems  of  interest  can  be  represented  (mathematically  modeled) 
linear  differential  equations.  For  exanple,  the  pitching  motion  of  an 
aircraft  in  flight  displays  motion  similar  to  a  mass-spring-danper  system  as 
shown  in  Figure  3.1. 


The  static  stability  of  the  aircraft  is  similar  to  the  spring,  the  moment  of 
inertia  about  the  y-axis  is  similar  to  the  mass,  and  the  airflow  (aerodynamic 
forces)  serves  to  danp  the  aircraft  motion.  Chapter  4  shows  that  stability 
derivatives  can  be  used  to  represent  the  static  stability  and  daiiping 
terms.  These  derivatives  are  C  and  C  .  In  this  chapter,  M,  K,  and  D  will 

ID  m  ** 

a  <3 

be  used  to  represent  mass,  spring,  and  danper  terms  respectively. 
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The  following  terms  will  be  used  extensively: 


Differential  Equation:  An  equation  relating  two  or  more  variables  in 
terms  of  derivatives. 

Independent  Variables:  Variables  that  are  not  dependent  on  other 

variables. 

Dependent  Vari^les:  Variables  that  are  dependent  on  other  variables. 
In  a  differential  equation,  the  dependent  variables  are  the  variables  on 
the  left-hand  side  of  the  equation  that  have  their  derivatives  taken  with 
respect  to  another  variable.  The  other  variable,  usually  time  in  our 
study,  is  the  independent  variable. 

Solution.  Any  function  without  derivatives  that  satisfies  a  differential 
equation. 

Ordinary  Differential  Equation.  A  differential  equation  with  only  one 
indepeiwent  variable. 

Partial  Differential  Equation.  A  differential  equation  with  more  than 
one  independent  variable. 

Order.  An  n*^**  derivative  is  a  derivative  of  order  n.  A  differential 
equation  has  the  order  of  its  highest  derivative. 

^gree.  The  expjnent  of  a  differential  term.  ihe  degree  of 
differential  equation  is  the  exponent  of  its  highest  order  derivative. 

Linear  Differential  Ekrjation.  A  differential  equation  in  vdiich  the 
dependent  variable  and  all  its  derivatives  are  only  first  degree,  and  the 
coefficients  are  either  constants  or  functions  of  the  independent 
variable. 

Linear  System.  Any  physical  system  that  can  be  described  \diich  satisfies 
a  differential  equation  of  order  n  vdiich  contains  n  arbitrary  constants. 

(^neral  Solution.  Any  function  without  derivatives  vAiich  satisfies  a 
differential  equation  of  order  n  which  contains  n  arbitrary  constants. 

3.2  BASIC  DIFFERENTIAL  EQUATION  SOLUTION 


Unfortunately,  there  is  no  general  method  to  solve  all  types  of 
differential  equations.  The  solving  of  a  differential  equation  involves 
finding  a  mathematical  expression  without  derivatives  v^ich  satisfies  the 
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differential  equation.  It  is  usually  much  easier  to  determine  vdiether  or  not 
a  candidate  solution  to  a  differential  equation  is  a  solution  than  to 
detertDine  a  likely  candidate.  For  exaitple,  given  the  linear  first  order 
differential  equation 


and  a  possible  candidate  solution 

y  =  ^  +  4x  +  C  (3.2) 

it  is  easy  to  differentiate  Equation  3.2  and  substitute  into  Equation  3.1  to 
see  if  Equation  3.2  is  a  solution  of  Equation  3.1.  The  derivative  of  Equation 
3.2  is 


x  +  4  (3.3) 

Substituting  Equation  3.3  into  Equation  3.1, 

(x  +  4)  -  X  =  4  (3.4) 

4  =  4 

Therefore,  Equation  3.2  is  a  solution  of  Equation  3.1. 

It  is  interesting  that,  in  general,  solutions  to  linear  differential 
equations  are  not  linear  functions.  Note  that  Equation  3.2  is  not  an  equation 
of  the  form 


y  =  mx  +  b  (3.5) 

\diich  represents  a  straight  line.  As  shown  in  Equation  3.2,  y  is  a  function 
of  X  and  x^ . 


iSiere  are  several  methods  in  use  to  solve  differential  equations.  Ih 
methods  to  be  discussed  in  this  chapter  are: 


1.  Direct  Integration 

2.  Separation  of  Variables 

3.  Exact  Differential  Integration 

4.  Integrating  Factor 

5.  Special  Procedures,  to  include  Operator  Techniques  and  Laplace 
Transforms. 

3.2.1  Direct  Integration 

Since  a  differential  equation  contains  derivatives,  it  is  sometimes 
possible  to  obtain  a  solution  by  anti-differentiation  or  integration.  This 
process  removes  the  derivatives  and  provides  arbitrary  constants  in  the 
solution.  For  exanple,  given 


rewriting 


integrating 


dy  -  xdx  =  4dx 


(3.1) 


(3.6) 


J  dy  -  J  xdx  =  J  4dx  +  C 


y  -  =  4x  +  C 


(3.7) 


or,  solving  for  y 


y  =  ^  +  4x  +  C 


(3.8) 


v4iere  C  is  an  arbitrary  constant  of  integration. 

Itofortunately,  application  of  the  direct  integration  process  fails  to 
work  in  many  cases. 
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3.2.2  Separation  of  Variables 


If  direct  integration  fails  for  a  first  order  differential  equation,  then 
the  next  step  is  to  try  to  separate  the  variables.  Direct  integration  may 
then  be  possible.  When  a  differential  equation  can  be  put  in  the  form 

f^(x)dx  +  f2(y)dy  =  0  (3.9) 

vdiere  one  term  contains  function  of  x  and  dx  only,  and  the  other  functions  of 
y  and  dy  only,  the  variables  are  said  to  be  separated.  A  solution  of  Equation 
3.9  can  then  be  obtained  by  direct  integration 

J  (x)dx  +  J  fj (y)dy  =  C  (3.10) 

where  C  is  an  arbitrary  constant.  Note,  that  for  a  differential  equation  of 
order  there  is  one  arbitrary  constant.  In  general,  the  number  of 
arbitrary  constants  is  equal  to  the  order  of  the  differential  equation. 

EXAMPLE 

+  3x  +  4 
ox  y  +  6 

(y  +  6)dy  =  (x^  +  3x  +  4)dx 

J  (y  +  6)dy  =  J  (x^  +  3x  +  4)dx  +  C 

^  +6y  =  |  +-^  +4x  +  C 

Not  all  first  order  equations  can  be  separated  in  this  fashion. 

3.2.3  Exact  Differential  Integration 

If  direct  integration,  or  direct  integration  after  separation  is  not 
possible,  then  it  still  may  be  possible  to  obtain  a  solution  if  the 
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differential  equation  is  an  exact  differential.  Associated  with  each  suitabl 
differentiable  function  of  two  variables  f(x,y),  there  is  an  expression  calle 
its  differential,  namely 


df  s 


that  can  be  written  as 


3x 


(3.11) 


df  =  M(x,y)dx  +  N(x,y)dy  =  0 

and  is  exact  if  and  only  if 

3y  3x 

If  the  differential  equation  is  exact,  then  for  all  values  of  C 

X  y 


(3.12) 


(3.13) 


M(x,y)dx  + 


N(x,y)cfy  =  C 


(3.14) 


is  a  solution  of  the  equation,  \diere  a  and  b  are  dumny  variables  o 
integration. 


EXAMPLE 

Show  that  the  equation 


(2x  +  3y  -2)dx  +  (3x  -  4y  +  l)dy  =  0  (3.15) 

is  exact  and  find  a  general  solution. 

Applying  the  test  in  Equation  3.13 

3M  3(2x  +  3y  -  2)  _  , 

3y  3y 

3N  _  3(3x  -  4y  +  1)  ^  , 

3x  3x  ^ 
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Since  the  two  partial  derivatives  are  equal,  the  equation  is  exact.  Its 
solution  can  be  found  by  means  of  Equation  3.14. 


X 


y 


J  (2x  +  3y  -  2)dx  + 
a 


J 

b 


(3x  -  4y  +  l)dy 


C 


The  integration  is  performed  assuming  y  is  a  constant  while  integrating  the 
first  term. 


X 


y 


(x^  +  3xy  -  2x) 


+  (3xy  -  2y^  +  y) 


C 


a 


b 


(x^  +  3xy  -  2x)  -  (a^  +  3ay  -  2a)  +  (3xy  -  2y^  +  y)  -  {3xb  -  2b^  +  b)  =  C 
x^  +  6xy  -  2x  -2y^  +  y  +  3ay  +  3xb  =  C  +  a^  -  2a  -  2b^  +  b  =  q  (3.16) 

The  same  result  can  be  obtained  with  less  algebra  and  probably  less 
chance  of  error  by  comparing  Equation  3.15  with  the  differential  form  in 
Equation  3.11. 


If 

dx 


dx  + 


- 


(3.11) 


(2x  +  3y  -  2)dx  +  (3x  -  4y  +  l)dy  =  0 


Conparing  these  two  equations, 

11=  2x  +  3y-2  =  0 

and 


||  =  3x  -  4y  +  1 


0 


(3.15) 


(3.17) 


(3.18) 
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Since  Equation  3.15  is  an  exact  differential,  then  Equations  3.17  and  3.18  ca^^ 
be  obtained  by  taking  partial  derivatives  of  the  same  function  f.  To  find  thl^r 
unknown  function  f,  first  integrate  Equations  3.17  and  3.18  assuming  that  y  is 
constant  vAien  integrating  with  respect  to  x  and  that  x  is  constant  v^en 
integrating  with  respect  to  y. 

f  =  x^  +  3xy  -  2x  +  f(y)  +  C  =  0  (3.19) 

f  =  3xy  -  2y^  +  y  +  f(x)  +  C  =  0  (3,20) 


Note  that  if  Equation  3.17  had  been  obtained  from  Equation  3.19,  e«iy  term  that 
was  a  function  of  y  only,  f(y),  and  any  constant  term,  C,  would  have 
disappeared.  Similarly,  obtaining  Equation  3.18  from  Equation  3.20,  the  f(x) 
and  C  terms  would  have  vanished.  By  a  direct  conparison  of  Equation  3.19  and 
3.20  the  total  fxjnction  f  can  be  determined. 


f  = 


x^  +  3xy  -  2x  -  2y^  +  y  +  C  =  0 


(3.21) 


Note  that  the  unknown  f(y)  term  in  Equation  3.19  is  (-2y^  +  y)  and  the  unknow^^ 
f(x)  term  in  Equation  3.20  is  x^  — 2x.  Redefining  the  constant  of  integration^^^^ 
Equation  3.21  can  be  written  as 


x^  +  3xy  -  2x  -  2y^  +  y  =  q  (3.16) 

and  was  obtained  earlier  integrating  using  dummy  variables  of  integration. 
3.2.4  Integrating  Factor 

When  none  of  the  above  procedures  or  techniques  work,  it  may  still  be 
possible  to  integrate  a  differential  equation  using  an  integrating  factor. 
When  some  unintegrable  differential  equation  is  multiplied  by  some  algebraic 
factor  which  permits  it  to  be  integrated  term  by  term,  then  the  algebraic 
factor  is  called  an  integrating  factor.  Determining  integrating  factors  for 
arbitrary  differential  equations  is  beyond  the  scope  of  this  course;  however, 
two  integrating  factors  will  be  introduced  in  later  sections  of  this  chapter 
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v^en  developing  operator  techniques  and  Laplace  transforms.  These  two  factors 
will  be  e"*  and  e"®^. 


3.3  FIRST  ORDER  EQUATIONS 

The  solution  to  a  first  order  linear  differential  equation  can  be 
obtained  by  direct  integration.  Consider  the  form 


g  +  R(x)y  =  0 

vdiere  R(x)  is  a  function  of  x  only  or  a  constant. 

^  +  R(x)dx  =  0 


(3.22) 

To  solve,  separate  variables 

(3.23) 


Integrating 

J  =  -  J  R(x)dx  +  C'  (3.24) 

vdiere 


Thus 


C'  =  In  C 

In  y  =  -  J  R(x)dx  +  In  C  (3.25) 


or 

y  =  Ce 

If  R(x)  is  a  constant,  R,  then 


y 


Ce" 


Rx 


(3.26) 

(3.27) 


From  this  result,  it  can  be  concluded  that  a  first  order  linear  differential 
equation  in  the  form  of  Equation  3.22  can  be  solved  by  simply  expressing  the 
solution  in  the  form  of  Equation  3.27. 


E3CAMPLE 


(3.28) 
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(3.29) 


then  the  solution  can  be  written  directly  as 


EXAMPLE 


is  in  the  form 


y  =  Ce" 


^  +  x^y  =  0 


^  +  R(x)y  =  0 


vdiich  has  the  solution 


y  =  Ce 


-  J  R(x)dx 


Therefore/  the  solution  to  Equation  3.30  can  be  obtained  directly 


x^dx 


y  =  Ce 


1 

■3  ^ 


y  =  Ce 


3.4  LINEAR  DIFFERENTIAL  EQUATIONS  AND  OPERATOR  TEf^TQTTKP 


(3.30) 


(3.22) 


(3.26) 


A  form  of  differential  equation  that  is  of  particular  interest 


"  .j.,®  j  n-l  +  •  •  •  + 


(3.31) 


If  the  coefficient  expression  ,  A^_j ,  .  .  .  »  Aq  are  all  functions  of  x 
only,  then  Equation  3.31  is  called  a  linear  differential  equation.  If  the 
coefficient  expressions  \  >  \  are  all  constants,  then  Equation  3.31 
is  called  a  linear  differential  equation  with  constant  coefficients. 
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EXAMPLE 


^2  d_y  ^2  ^  +  xy  =  sin  x 
dx^ 

is  a  linear  differential  equation. 

EXAMPLE 


=  e^ 

dx^  ^ 


is  a  linear  differential  equation  with  constant  coefficients.  Linear 
differential  equations  with  constant  coefficients  occur  frequently  in  the 
analysis  of  physical  systems.  Mathematicians  and  enqineers  have  developed 
single  and  effective  techniques  to  solve  this  type  of  equation  by  using  either 
"classical"  or  operational  methods.  when  attenpting  to  solve  a  linear 
differential  equation  of  the  form 


dx" 


.  + 


AqY  = 


(3.32) 


it  is  helpful  to  first  examine  the  equation 


(3.33) 


Equation  3.33  is  the  same  as  Equation  3.32  with  the  right-hand  side  set  equal 
to  zero.  Equation  3.32  is  known  as  the  general  equation  and  Equation  3.33  as 
the  conplementary  or  homogeneous  equation.  Solutions  of  Equation  3.33  possess 
a  useful  property  known  as  superposition,  vdiich  may  be  briefly  stated  as 
follows;  Suppose  y^ (x)  and  y2(x)  are  distinct  solutions  of  Equation  3.33. 
Then  any  linear  combination  of  y^ (x)  and  y2(x)  is  also  a  solution  of  Equation 
3.33.  A  linear  combination  would  be  C^y^(x)  +  C2y2(x). 
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EXAMPLE 


—  0 

It  can  be  verified  that  y^Cx)  =  e^*  is  a  solution,  and  that  Yj  (x)  =  e^*  is 

another  solution  vdiich  is  distinct  from  yj(x).  Using  superposition,  then, 
y(x)  =  +  Cje^*  is  also  a  solution. 

Equation  3.32  may  be  interpreted  as  representing  a  physical  system  where 
the  left  side  of  the  equation  describes  the  natural  or  designed  state  of  the 
system,  and  vdiere  the  right  side  of  the  equation  represents  the  input  or 
forcing  function. 

Ihe  following  line  of  reasoning  is  used  to  find  a  solution  to  Equation 

3.32; 


1.  A  general  solution  of  Equation  3.32  must  contain  n  arbitrary 
constants  and  must  satisfy  the  equation. 

2.  The  following  statements  are  justified  by  experience; 

a.  It  is  reasonably  straightforward  to  find  a  solution  to  the 
complementary  Equation  3.33,  containing  n  arbitrary 
constants.  Such  a  solution  will  be  called  the  transient 
solution.  Physically,  it  represents  the  response  present 
in  the  system  regardless  of  input. 

b.  There  are  varied  techniques  for  finding  the  solution  of  a 
differential  equation  due  to  a  forcing  function.  Such 
solutions  do  not,  in  general,  contain  arbitrary  constants. 
This  solution  will  bs  called  the  particular  or  steady 
state  solution. 

3.  If  the  transient  solution  vdiich  describes  the  response  already 

existing  in  the  system  is  added  to  the  response  due  to  the 

forcing  fmction,  it  would  appear  that  a  solution  so  written 
would  blend  the  two  responses  and  describe  the  total  response 
of  the  system  represented  by  Equation  3.32.  In  fact,  the 
definition  of  a  general  solution  is  satisfied  under  such  an 
arrangement.  This  is  sinply  an  extension  of  the  principle  of 
superposition.  The  transient  solution  contains  the  correct 

number  of  arbitrary  constants,  and  the  particular  solution 
guarantees  that  the  combined  solutions  satisfy  the  general 

Equation  3.32.  A  general  solution  of  Equation  3.32  is  then 
given  by 
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y 


(3.34) 


=  yt  +  yp 

vrtiere  is  the  transient  solution  and  is  the  particular 
solution. 

3.4.1  Transient  Solution 

Equation  3.28  is  a  conplementary  or  homogeneous  first  order  linear 
differential  equation  with  constant  coefficients.  A  quick  and  sinple  method 
of  solving  this  equation  was  found.  The  solution  was  always  of  exponential 
form;  hopefully,  solutions  of  higher  order  equations  of  the  same  family  take 
the  same  form. 


^  +  2y  =  0  (3.28) 

Next,  a  second  order  differential  equation  with  constant  coefficients 
will  be  examined  to  determine  if  the  candidate  solution 

y  =  e"*  (3.35) 

is  a  solution  of  the  equation 

ay"  +  by'  +  cy  =  0  (3.36) 

v^en  the  prime  notation  indicates  derivatives  with  respect  to  x.  That  is, 

d^y 

dx^ 

Substituting 

y  =  e"" 

am^e"*  +  bme”*  +  ce"**  =  C 

or 

(am^  +  bm  +  c)e'"*  =  0. 

Since 

e’”'  ^  0 

am^  +  tan  +  c  =  0 


(3.37) 

(3.38) 
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(3.39) 


and,  using  the  quadratic  formula 


®“i,2  =  -b  ±  Vb^  -  4ac 


Substituting  these  values  into  the  assumed  candidate 
when  m^  and  are  defined  by  Equation  3.40. 

m^x  nijX 

Yt  ®  ^1®  +  CjO 


(3.40) 


solution,  it  is  a  solution 


(3.41) 


Equation  3.41  represents  a  transient  solution  since  there  is  no  forcing 
function  in  Equation  3.36.  When  working  numerical  problems,  it  is  not 
necessary  to  take  the  derivatives  of  e"* .  This  will  be  true  for  any  order 
differential  equation  with  constant  coefficients.  From  the  foregoing,  it  is 
seen  that  the  method  for  first  order  conplementary  equations  has  been  extended 
to  higher  order  conplementary  or  homogeneous  equations.  Again  an  integration 
problem  has  been  traded  for  an  algebra  problem  (solving  Equation  3.39  for  m's). 

There  are  four  possibilities  for  m^  and  m^ ,  and  each  is  discussed  below. 

3. 4. 1.1  Case  1;  Roots  Real  and  Unequal.  If  m^  and  nij  are  real  and  unequal, 
the  desired  form  of  solution  is  just  as  given  Equation  3.41. 

EXAMPLE  ^ 

Given  the  homogeneous  equation 


— ^  +  4  -  12y  =  0, 

dx^ 


rewriting  in  operator  form  vAiere 


(m^  +  4in  -  12)y  =  0. 


Solving  for  the  values  of  m, 


gives 


m  +  4m  -  12  =  0 
-4  +  Vl6  +  48 


-4  +  8 


m  = 


=  -6,2 


and  the  required  transient  solution  is 


Yt  =  +  Cje^V 


3. 4. 1.2.  Case  2:  Roots  Real  and  Equal.  If  m^  and  nij  are  real  and  equal,  an 
alternate  form  of  solution  is  required. 


EXAMPLE 

Given  the  homogeneous  equation 


£y  _ 

dx^ 


4  ^  +  4y  =  0, 


(3.42) 


rewriting  in  operator  form 


(m  -  4m  +  4)y  =  0. 

Solving  for  the  values  of  m. 


4  +  Vl6  -  16 


m  = 


=  -  =  2 
2  ^ 


or  m  =  2.  But  this  gives  only  one  value  of  m,  and  two  values  of  m  are 
required  to  result  in  a  solution  of  the  form  of  Equation  3.41  v4iich  has  two 
arbitrary  constants.  The  operator  expression 


can  also  be  written 


m  -  4m  +  4  =  0 


(m  -  2f  =0 


3.15 


or 


(m  -  2)  (m  -  2) 


now  a  repeated  polynomial  factor  resulting  in  two  (repeated)  roots 

m  =  2,2. 


Writing  the  solution  in  the  form  of  Eguation  3.41  \dien  the  roots  are  repeated 
does  not  give  a  solution  because  the  two  arbitrary  constants  can  be  combined 
into  a  single  arbitrary  constant  as  shown  below. 


Yt  =  +  c^e^*  =  (c^  +  C2)e^’‘  =  CjO^* 

To  solve  this  problem  one  of  the  arbitrary  constants  is  multiplied  by  x.  The 
solution  now  contains  two  arbitrary  constants  vhich  cannot  be  combined,  and  it 
is  easily  verified  that 


y^  =  +  Cj*® 

is  a  transient  solution  of  Equation  3.42. 

3. 4. 1.3  Case  3;  Roots  Purely  Imaginary. 


EXAMPLE 

Given  the  homogeneous  equation 

32  _ 


dx^ 


y  =  0, 


rewriting  in  operator  form 


Solving, 


(m^  +  l)y  =  0. 

0  +  Vo  -  4 


m  = 


- 


3.16 


In  most  engineering  work  1  is  given  the  symbol  j .  ( In  mathematical  texts 

it  is  denoted  i . ) .  Now, 

m  =  +  j 

and  the  solution  is  written 

=  c,e^"  +  c^e"^*  (3.43) 

This  is  a  perfectly  good  solution  from  a  mathematical  standpoint,  but  Euler's 
identity  can  be  used  to  put  the  solution  in  a  more  useable  form. 

=  cos  X  +  jsin  x  (3.44) 

This  equation  can  be  restated  in  many  ways  geometrically  and  analytically,  and 
can  be  verified  by  adding  the  series  expansion  of  cos  x  to  the  series 
expansion  of  j  sin  x .  Now  Equation  3.43  may  be  expressed 

y^  =  Cj(cos  X  +  jsin  x)  +  Cjtcos  (-x)  +  jsin(-x)] 

“  (Ci  +  Cj )cos  X  +  j(Cj  -  Cj )sin  X  (3.45) 

or  without  loss  of  generality 

y^  =  C3  cos  X  +  C4  sin  x  (3.46) 

An  equivalent  expression  to  Equation  3.46  is 


cj  +  c/  ( 


c  ^  +  c  ^ 


cos  X  + 


sin  x) 


(3.47) 


c,  and  c. 

3  4 


If  the  arbitrary  constants 


are  related  as  shown  in  Figure  3.2, 


then 


FIGURE  3.2.  DEFINITIC»I  OF  C,  AND  c. 

3  4 


and 


sin  <f> 


cos  <|» 


=  * 

vdiere  A  and  are  also  arbitrary  constants,  Equation  3.47  becomes 

=  A  (sin  <j»  cos  x  +  cos  4>  sin  x) 
or  using  a  common  trigononetric  identity 

y^  =  A  sin(x  +  ♦) 

Note  also  that  Equation  3.48  could  be  written  in  the  equivalent  form 

y,.  =  A  cos(x  -  0) 


vdiere 


e  =  90"  -  <f» 


(3.48) 


(3.49) 
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TO  svunmarize,  if  the  roots  of  the  operator  polynomial  are  purely  imaginary, 
they  will  be  numerically  equal  but  opposite  in  sign,  and  the  solution  will 
have  the  form  of  Equation  3,46,  3.48,  or  3.49. 

EXAMPLE 

Given  the  homogeneous  equation 

ly  .  0 

dx^ 

rewriting  in  operator  form 

(m^  +  4)y  =  0 


vdiich  gives  the  roots 

m  =  +  2j 


Alternate  solutions  can  immediately  be  written  as 

y.  =  c-  cos  2x  +  c  sin  2x 
or 

y,.  =  A  sin(2x  +  <()) 

vdiere  Cj ,  c^ ,  A,  and  are  arbitrary  constants. 

3. 4. 1.4  Case  4:  Roots  Conplex  Conjugates. 

EXAMPLE 


Given  the  homogeneous  equation 

^  +  2  ^  +  2y  =  0 

dx^  ^ 

rewriting  in  operator  form 

(m^  +  2m  +  2)y  =  0 

Solving  gives  a  coitplex  pair  of  roots 


m  = 


-2  ±  \/4  -  8 

2 


-1  ±)fT 
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or 


m  =  -1  +  j,  -1  -  j 
The  solution  can  be  written 


y  —  c  e^~^  ^  ~  j)* 


Factoring  out  the  exponential  term  gives 


y,  =  e-"(c,e^"  +  c^e"^") 


ot f  using  the  results  from  Equations  3.46  and  3.48,  alternate  solutions  can 
inonediately  be  written  as 


yt 

=  COS  X  +  c^  sin  x) 

(3.50) 

Yt 

=  e“*  A  sin(x  +  4>) 

(3.51) 

3.4.2  Particular  Solution 

^e  particular  solution  to  a  linear  differential  equation  can  be  obtained 
by  the  method  of  undetermined  coefficients.  This  method  consists  of  assuming 
a  solution  of  the  same  general  form  as  the  input  (forcing  function),  but  with 
undetermined  constant  coefficients.  Substitution  of  this  assumed  solution 
into  the  differential  equation  enables  the  coefficients  to  be  evaluated.  The 
method  of  undetermined  coefficients  applies  when  the  forcing  function  or  input 
is  a  polynomial,  or  of  the  form 


sin  ax,  cos  ax,  e“* 


or  combinations  of  sums  and  products  of  these  terms.  The  general  solution  to 
the  differential  equation  with  constant  coefficients  is  then  given  by  Equation 
3.34, 

y  =  Yt  +  Yp  (3.34) 
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vAiich  is  the  summation  of  the  solution  to  the  complementary  ecjuation 
(transient  solution),  plus  the  particular  solution. 

Consider  the  equation 


a  — f  ^  ^  +  cy  =  f(x)  (3.52) 

dx 

The  particular  solution  vrtiich  results  from  a  given  input,  f(x),  can  be  solved 
for  using  the  method  of  undetermined  coefficients.  The  method  is  best 
illustrated  by  considering  examples. 

3. 4. 2.1  Constant  Forcing  Functions. 

EXAMPLE 


dx^ 


+  3y  =  6 


(3.53) 


The  input  is  a  constant  (trivial  polynomial),  so  a  solution  of  form  y  =  K 

, 

IS  assumed. 

Then 


dK 


dx  dx 


=  0 


and 


d^y. 


dx^ 


i  =  0 

dx^ 


Substituting  into  Equation  3.53, 

0  +  4(0)  +  3K  =  6 


yp  =  K  =  2 
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Therefore,  =  2  is  a  particular  solution, 

solved  using  operator  form 


The  homogeneous  equation  can 


or 


— ^  +  4  ^  +  3y  =  0 
dx^ 

(m^  +  4m  +  3)y  =  0 
m  =  -1,  -3 


and  the  transient  solution  can  be  written  as 


(3.54) 


yt 


c^e“*  + 


(3.55) 


The  general  solution  of  Equation  3.53  is 


y 


c^e  *  +  c^e”^* 

transient 

solution 


+  2 
^ - ^ 

particxilar 

(or  steady  state) 

solution 


3. 4. 2. 2  Polynomial  Forcing  Function. 
EXAMPLE 


(3.56) 


+  4  ^  +  3y  =  +  2x 

dx^  ^ 


(3.57) 


The  form  of  f(x)  for  Equation  3.57  is  a  polynomial  of  second  degree,  so  a 
particular  solution  for  y^  of  second  degree  is  assumed: 

y^  =  Ax^  +  Bx  +  C 


Then 


=  2Ax  +  B 
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and 


dx^ 

Substituting  into  Equation  3.57, 

(2A)  +  4(2Ax  +  B)  +  3{Ax^  +  Bx  +  C)  =  x^  +  2x 
or 

(3A)x^  +  (8A  +  3B)x  +  (2A  +  4B  +  3C)  =  x^  +  2x 

Equating  like  powers  of  x, 

xS  3A  =  1 

A  *  1/3 

x;  8A  +  3B  =  2 

3B  =  2  -  8/3 

B  =  -  2/9 

X® :  2A  +  4B  +  3C  =  0 

3C  =  8/9  -  2/3 
C  =  2/27 

Therefore, 

Yp  =  1/3  x^  -  2/9  X  +  2/27 
The  total  general  solution  of  Equation  3.57  is  given  by 

y  =  c^e"*  +  c^e'^*  +  1/3  x^  -  2/9  x  +  2/27  (3.58) 

since  the  transient  solution  is  Equation  3.55.  As  a  general  rule,  if  the 
forcing  function  is  a  polynomial  of  degree  n,  assume  a  polynomial  solution  of 
degree  n. 
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3. 4. 2. 3  EXPCaJENTIAL  FORCING  FUNCTION. 
EXAMPLE 


^  ^  (3.59 

dx 

The  forcing  frinction  is  so  assume  a  solution  of  the  form 

Yp  =  Ae^* 

^  (Ae^*)  =  2Ae^* 

-'(Ae^M  =  4Ae"* 

dx^ 

Substituting  into  Equation  3.59, 

4Ae^’‘  +  4(2Ae^’‘)  +  3(Ae^*)  =  e^’' 

e^*(4A  +  8A  +  3A)  =  e^*  ( 

The  coefficients  on  both  sides  of  the  equation  must  be  the  same.  Therefore, 
4A+8A+3A  =  1,  or  15a  =  1,  and  A  =  1/15.  The  particular  solution  of 
Equation  3.59  then  is  =  1/15  e^*.  The  transient  solution  is  still 

Equation  3.55.  A  final  exanple  will  illustrate  a  pitfall  sometimes 
encountered  using  this  method. 
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3. 4.2.4 
EXAMPLE 


^  +  4  ^  +  3y  =  e“* 
dx^  ^ 

The  forcing  function  is  e"’‘,  so  assume  a  solution  of  the  form  v 
Then 

^  (ite- )  =  -Ae- 

and 

—  (Ae“’‘)  =  Ae"* 

dx^ 

Substituting  into  Equation  3.60, 

Ae"*  +  4(-Ae“* )  +  3(Ae~* )  =  e"* 

(A  -  4A  +  3A)e"*  =  e"* 

(0)e"’‘  =  e"’^ 

Obviously,  this  is  an  incorrect  statement.  To  locate  the  difficulty,  the 
procedure  to  solve  differential  equations  will  be  reviewed. 

To  solve  an  ecfuation  of  the  form 

(m  +  a)  (m  +  b)y  =  e”* 

solve  the  homogeneous  equation  to  get 

(m  +  a)(m  +  b)y  =  0 

m  =  -a,  -b 

yt  =  Cie  +  c^e 


3.25 


is  assumed  for  a  particvilar  solution,  then 


If  y  =  Ae" 


y  =  Yt  +  yp  “  C^e  “*  +  +  Ae"“’‘  =  (c^  +  A)e"*’‘  +  c^e 


_  “bx 


-  C3e-“*  + 


Yt 


However,  is  the  solution  only  vdien  the  right  side  of  the  equation  is  zero, 
and  will  not  solve  the  equation  vdien  there  is  a  forcing  function  of  the  form 
given.  Asstmiing  a  particular  solution  of  the  form 

y„  =  Axe"** 

will  lead  to  a  solution,  then 

y  =  Yp  +  Yt  =  c^e"**  +  Cje"”*  +  Axe"“*  =  (c^  +  Ax)e"**  +  c^e"*’*  y^ 

Similarly,  the  equation 

(m  +  aj)(m  -  aj)y  =  sin  ax  ^ 

has  the  transient  solution 

y^  =  Cj  sin  ax  +  Cj  cos  ax 

If  Yp  =  A  sin  ax  +  B  cos  ax  is  assximed  for  a  particular  solution,  then 

y  =  Yt  +  Yp  =  +  A)sin  ax  +  (Cj  +  b)cos  ax 

y  =  C3  sin  ax  +  c^  cos  ax  =  y^. 

vdiich,  as  in  the  previous  exanple,  does  not  provide  a  solution  ^en  there  is  a 
forcing  function  of  the  form  given.  But,  assuming  a  solution  of  the  form 


yp  =  Ax  sin  ax  +  Bx  cos  ax 
does  lead  to  a  solution 


y 


(c^  +  Ax)sin  ax  +  (c^  +  Bx)cos  ax 


Yt 
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Continuing  with  the  solution  of  Equation  3.60,  a  valid  solution  can  be  foiand 
by  assuming  =  Axe~*  ,  then 

^  (Axe~*  )  =  A(-xe“*  +  e"*  ) 

and 

j2 

— ,  (Axe"* )  =  A(xe"*  -  2e"*  ) 

dx^ 

Substituting  into  Equation  3.60, 

A(xe"*  -  2e"*)  +  4A(-xe"*  +  e"*)  +  3(Axe"*)  =  e"* 

(A  -  4A  +  3A)xe"*  +  (-2A  +  4A)e"*  =  e“* 

(0)xe"*  +  2Ae"*  =  e"* 

and 


A  =  1/2 

Thus, 

yp  =  (l/2)xe-* 

is  a  particular  solution  of  Equation  3.60,  and  the  general  solution  is  given 


y  =  c^e"*  +  c^e"^*  +  1/2  xe"* 

The  key  to  successful  application  of  the  method  of  undetermined  coefficients 
is  to  assxame  the  proper  form  for  a  trial  or  candidate  particular  solution. 
Table  3.1  summarizes  the  results  of  this  discussion.  When  f(x)  in  Table  3.1 
consists  of  a  sum  of  several  terms,  the  appropriate  choice  for  y  is  the  sum 
of  y^  expressions  corresponding  to  these  terms  individually.  Whenever  a  term 
in  any  of  the  Yp'^  listed  in  Table  3.1  duplicates  a  term  already  in  the 
conplementary  function,  all  terms  in  that  y^  must  be  multiplied  by  the  lowest 
positive  integral  power  of  x  sufficient  to  eliminate  the  duplication. 
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3.4.3  Solving  For  Constants  of  Integration 

As  discussed  previously,  the  number  of  arbitrary  constants  in  the 
solution  of  a  linear  differential  equation  is  equal  to  the  order  of  the 
equation.  The  constants  of  integration  can  be  determined  by  initial  or 
boundary  conditions.  That  is,  to  solve  for  the  constants  the  physical  state 
(position,  velocity,  etc.)  of  the  system  must  be  known  at  some  time.  The 
number  of  initial  or  boxmdary  conditions  given  must  equal  the  number  of 
constants  to  be  solved  for.  Many  times  these  conditions  are  given  at  time 
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equal  to  zero,  in  vdiich  case  they  are  called  initial  conditions.  A  system 
vdiich  has  zero  initial  condition,  i.e.,  initial  position,  velocity,  and 
acceleration  all  equal  to  zero,  is  frequently  called  a  quiescent  system. 

The  arbitrary  constants  of  the  solution  must  be  evalxiated  from  the  total 
general  solution,  that  is,  the  transient  plus  the  steady  state  solution.  The 
method  of  evaluating  the  constants  of  integration  will  be  illustrated  with  an 
example . 

EXAMPLE 


*x*  +  4x  +  13x  =  3  (3.61) 

where  the  dot  notation  indicates  derivatives  with  respect  to  time,  that  is,  x 
=  dx/dt,  X  =  d  x/dt  .  The  initial  conditions  given  are  x(0)  =  5,  and 
x(0)  =  8.  The  transient  solution  is  given  by 

m^  +  4m  +  13  =  0 

m  =  -2  +  \/4  -  13  =  -2  +  j3 

x^.  =  e~^*^(A  cos  3t  +  B  sin  3t) 

Assimte  the  particular  solution  of  the  form 

Xp  =  D 


Substituting  into  Equation  3.61,  D  =  3/13  for  the  total  general  solution 
x(t)  =  e“^^(A  cos  3t  +  B  sin  3t)  +  3/13 
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To  solve  for  A  and  B,  the  initial  conditions  specified  above  are  used. 

x(0)  =  5  =  A  +  3/13 
or 

A  =  62/13 

Differentiating  the  total  general  solution, 

x(t)  =  e“^*'[(3B  cos  3t  -  3A  sin  3t)  -  2e“^^{A  cos  3t  +  B  sin  3t)] 
Substituting  the  second  initial  condition 

x(0)  =  8  *  3B  -  2A 


B  = 


76 

13 


Therefore,  the  conplete  solution  to  Equation  3.61  with  the  given  initial 
conditions  is 


x{t)  =  e'^’^E (62/13)  cos  3t  +  (76/13)  sin  3t]  +  3/13 

First  and  second  order  differential  equations  have  been  discussed  in  soim 
detail.  It  is  of  great  in^rtance  to  note  that  mai^  higher  order  system? 
quite  naturally  deconpose  into  first  and  second  order  systems.  For  example, 
the  study  of  a  third  order  equation  (or  system)  may  be  conducted  examining 
®  first  and  a  second  order  system,  a  fourth  order  system  analyzed  by  examining 
two  second  order  systems,  etc.  All  these  cases  are  handled  solving  the 
characteristic  equation  to  get  a  transient  solution  and  then  obtaining  the 
particular  solution  by  any  convenient  method. 

A  few  remarks  are  appropriate  regarding  the  second  order  linear 
differential  equation  with  constant  coefficients.  Although  the  equation  is 
interesting  in  its  own  right,  it  is  of  particular  value  because  it  is  a 
mathematical  model  for  several  problems  of  physical  interest. 
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a  ^  +  cy  =  f(x)  niathematical  mcxJel 

dx^  ^ 

X  dx 

M  —  ^  3t  describes  a  mass  spring 

danper  system 

L  — ^  ®  ^  C  ”  describes  a  series  LRC 

electrical  circuit 


(3.62) 

(3.63) 


(3.64) 


Equations  3.62,  3.63,  and  3.64  are  all  the  same  mathematically,  but  are 
expressed  in  different  notation.  Different  notations  or  symbols  are  enployed 
to  emphasize  the  physical  parameters  involved,  or  to  force  the  solution  to 
appear  in  a  form  that  is  easy  to  interpret.  In  fact,  the  similarity  of  these 
last  two  equations  may  suggest  how  one  might  design  an  electrical  circuit  to 
simulate  the  operation  of  a  mechanical  system. 


3.5  APPLICATIONS  AND  STftNDftRD  FORMS 

to  this  point,  differential  equations  in  general  and  linear 
differential  equations  with  constant  coefficients  have  been  considered. 
Methods  for  solving  first  and  second  order  equations  of  the  following  type 
have  been  developed: 


a  +  bx  =  f(t)  (3.65) 

a^  +  b|^  +  cx  =  f(t)  (3.66) 

dt 

These  two  equations  are  mathematical  models  or  forms.  These  same  forms  may  be 
used  to  describe  diverse  physical  systems.  This  section  will  concentrate  on 
the  transient  response  of  the  systems  under  investigation. 


3 


(3.67) 


3.5.1  First  Order  Equation 
EXAMPLE 


4x  +  X  = 

Physically,  x  can  represent  distance  or  displacement,  vdiere  t  is  used  to 
represent  time.  The  transient  solution  can  be  found  from  the  homogeneous 
equation. 


4x  +  X  =  0 
(4m  +  l)x  =  0 
4m  +  1  =  0 
m  =  -1/4 

Hius 

x^  =  ce-"/^ 

The  particular  solution  is  foiond  by  assuming 


Substitute 


or 


^P 

ar  = 


0 


A  =  3 
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The  total  general  solution  is  then 


X  =  +  3  (3.68) 

The  first  term  on  the  right  of  Equation  3.68  represents  the  transient 
response  of  the  physical  system  described  by  Equation  3.67,  and  the  second 
term  represents  the  steady  state  response  if  the  transient  decays.  A  term 
useful  in  describing  the  physical  effect  of  a  negative  exponential  term  is 
time  constant  \diich  is  denoted  by  t.  The  time  constant  is  defined 

T  =  -  - 

m 

Thus,  Equation  3.68  could  be  rewritten  as 


X  =  +  3  (3.69) 


\diere  t  =  4. 


Note  the  following  points: 

1.  The  time  constant  is  discussed  only  if  m  is  negative.  If  ro  is 
positive,  the  exponent  of  e  is  positive,  and  the  transient 
solution  will  not  decay. 

2.  If  m  is  negative,  t  is  positive. 

3.  X  is  the  negative  reciprocal  of  m,  so  that  small  numerical 
values  of  m  give  large  numerical  values  of  t  (and  vice  versa). 

4.  The  value  of  x  is  the  time,  in  seconds,  required  for  the 
displacei^nt  to  decay  to  1/e  of  its  original  displacement  from 
equilibrium  or  steady  value.  To  get  a  better  understanding  of 
this  statement,  examine  Equation  3.69 


Thus,  vdien  t  -  x,  the  exponential  portion  of  the  solution  has  decayed  to  1/e 
of  its  original  displacement  as  shown  in  Figure  3.3. 
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X 


FIGURE  3.3,  EXAMPLE  OF  FIRST  ORDER  EXPONENTIAL 
DECAY  WITH  AN  ARBITRARY  OXJSTANT 


Other  measures  of  time  are  sometimes  used  to  describe  the  decay  of  the 
exponential  of  a  solution.  If  T^  is  used  to  denote  the  time  it  takes  for  the 
transient  to  decay  to  one-half  its  original  amplitude,  then 

T^  *  0.693  T  (3.71) 

This  relationship  can  be  easily  shown  by  investigating 

X  =  c^e"“^  +  Cj  (3.72) 

By  definition,  t  =  1/a.  T^  is  the  value  of  t  at  vrtiich  x^  =  1/2  x^(0). 
Solving 


=  c^e 


-at 


-aT 

1/2  x^  ( 0 )  =  1/2  Cj  =  Cj  e  ^ 


=  1/2 


-In  1/2  =  aT^ 
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0.693T 


T  =  -In  1/2  ^  0.693  , 

1  a  a  ■ 

The  solution  of  Equation  3.67  can  be  completed  by  specifying  a  boiindary 
condition  and  evaluating  the  arbitrary  constant.  Let  x  =  0  at  t  =  0. 

X  =  ce“'^^^  +  3 
x(0)  =  0  =  c  +  3 

c  =  -3 

The  complete  solution  for  this  boundary  condition  is 

/ 

X  *  +  3 

as  shown  in  Figure  3.4. 


X 


FIGURE  3.4.  EXAMPLE  OF  FIRST  ORDER  EXPONENTIAL  DECAY 
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3.5.2  Second  Order  Equations 

Consider  an  equation  of  the  form  of  Equation  3.66 

+CX  -  f(t) 

dt“  “ 

the  characteristic  equation  can  be  written  in  operator 

am^  +  bm+c  =  0  (3.39) 

vrtiere  roots  can  be  represented  by 

-b  +  V  b^  -  4ac 

=  - 2i -  (3.40) 

These  quadratic  roots  determine  the  form  of  the  transient  solution.  The 
physical  in^lications  of  solutions  for  various  values  of  m  will  now  be 
discussed. 


As  discussed  earlier, 
notation  as 


3. 5. 2.1  Case  1 ;  Roots  Real  and  Unequal . 

unequal,  the  transient  solution  has  the  form 


When  the  roots  are  real  an( 


(3.73) 


When  m^  and  m^  are  both  negative,  the  system  decays  and  there  will  be  a 
time  constant  associated  with  each  exponential  as  shown  in  Figure  3.5. 


x,=  c,e  I'+Cje;  2‘ 
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FIGURE  3.5.  SECOND  ORDER  TRANSIENT  RESPONSE 

WITH  REAL,  UNEQUAL,  NEGATIVE  ROOTS 

When  nij^  or  nij  (or  both)  is  positive,  the  system  will  generally  diverge  as 
shown  in  Figures  3.6  and  3.7. 
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FIGURE  3.6.  SECCWD  ORDER  TRANSIENT  FIGURE  3.7.  SECOND  ORDER  TRANSIENT 
RESPC»ISE  WITH  ONE  RESPONSE  WITH  REAL, 

POSITIVE  AND  ONE  NEGATIVE  UNEQUAL,  POSITIVE  ROOTS 

REAL,  UNEQUAL  ROOTS 


3. 5. 2. 2  Case  2;  Roots  Real  and  Equal.  Nhen 
solution  has  the  fonn 


= 


the 


transie 


Xt  *  c^e"’^  +  Cjte"^  (3.74) 

When  m  is  negative,  the  system  will  usually  decay  as  shown  in  Figure  3.8. 
If  m  is  very  small,  the  system  may  initially  exhibit  divergence. 


X 


FIGURE  3.8.  SECCM)  ORDER  TRANSIENT  RESPONSE 
WITH  REAL,  EQUAL,  NEGATIVE  ROOTS 


When  m  is  positive,  the  system  will  diverge  much  the  same  way  as  shown  in 
Figure  3.7. 


3. 5. 2. 3  Case  3;  Roots  Purely  Imaginary.  When  m  =  +  jk,  the  transient 

solution  has  the  form 


or 


or 


x^,  =  Cj  sin  kt  +  Cj  cos  kt 

x^  =  A  sin(kt  +  +) 
x^  ®  A  cos(kt  +  0) 


(3.75) 

(3.76) 

(3.77) 
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I 


The  system  executes  oscillations  of  constant  amplitude  with  a  frequency  k  as 
shown  in  Figure  3.9. 


X 


FIGURE  3.9.  SECCM)  ORDER  TRANSIENT  RESPONSE 
WITH  IMAGINARY  ROOTS 


3. 5. 2. 4  Case  4;  Roots  Complex  Conjugates.  When  the  roots  are  given  by 
m  =  k^  +  jkj ,  the  form  of  the  transient  solution  is 


or 


or 


k  t 

el  (Cj  cos  kjt  +  Cj  sin  kjt) 

(3.78) 

Ae^i^  sinCkjt  +  <f>) 

(3.79) 

Ae^i^  cos(k2t  +  0) 

(3.80) 

The  system  executes  periodic  oscillations  contained  in  an  envelope  given 


Iy  X  =  + 

When  kj  is  negative,  the  system  decays  or  converges  as  shown  in  Figure 
3.10.  When  k^^  is  positive,  the  system  diverges  as  shown  in  Figure  3.11. 
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FIGURE  3.10.  SECXM)  ORDER  CX3NVERGENT  FIGURE  3.11.  SECCM)  ORDER  DIVERGENT 
TRANSIENT  RESPONSE  WITH  TRANSIENT  RESPC»JSE  WITH 

CXOTPLEX  CONJUGATE  ROOTS  COMPLEX  CCKJUGATE  ROOTS 


The  discussion  of  transient  solutions  above  reveals  only  part  of  the  pictura 
presented  by  Equation  3.66.  The  input  or  forcing  function  is  still  left 
consider,  that  is,  f(t).  In  practice,  a  linear  system  that  possesses  a 

divergence  (without  input)  may  be  changed  to  a  damped  system  by  carefully 
selecting  or  controlling  the  input.  Conversely,  a  nondivergent  linear  system 
with  weak  danping  may  be  made  divergent  certain  types  of  inputs.  Chapter 
13,  Linear  Control  Theory,  will  examine  these  problems  in  detail. 


3.5.3  Second  Order  Linear  Systems 

Consider  the  physical  model  shown  in  Figure  3.12.  The  system  consists  of 
an  object  suspended  by  a  spring,  with  a  spring  constant  of  K.  The  mass 
represented  by  M  may  move  vertically  and  is  subject  to  gravity,  input,  and 
damping,  with  the  total  viscous  damping  constant  equal  to  D. 
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FIGURE  3.12.  SECCM)  ORDER  MASS,  SPRING, 
DAMPER  SYSTEM 


The  equation  for  this  system  is  given  by 

+  Dx  +  Kx  =  f(t) 

The  characteristic  equation  in  operator  notation  is  given  by 

Mm^  +  Dm  +  K  =  0 

The  roots  of  this  equation  can  be  written 


(3.81) 

(3.82) 

(3.83) 

(3.84) 
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For  sinplicity,  eind  for 
constants  are  defined 


reasons  that  will 


be  obvious  later  thr 


C 


(3.85) 


the  term  C  is  called  the  damping  ratio,  and  is  a  value  v^ich  indicates  the 
danping  strength  in  the  system. 


(0 


n 


(3.86) 


w  is  the  undanped  natural  frequency  of  the  system.  This  is  the  frequency  at 
which  the  system  would  oscillate  if  there  were  no  danping  present. 

“d  =  y  1  -  ^  (3.87) 


is  the  damped  frequency  of  the  system,  it  is  the  frequency  at  which  the 
system  oscillates  v^en  a  damping  ratio  of  C  is  present. 

Substituting  the  definitions  of  C  and  into  Equation  3.84  gives 


^  ±  :j"n 


(3.88)i 


With  these  roots,  the  transient  solution  becomes 


m  t  m  t 
=  c^e  +  c^e  ^ 


(3.89) 


vdiich  can  be  written  as 


-Cw  t 


X.  =  e 


"  (c,  cos  0)  Vl  -  t  +  c.  sin  w  V  1  - 


^  t) 


(3.90) 


or 


-Cco  t  y - 

X  =  A  e  "  sin  (w  V 1  -  t  +  <)>) 


(3.91) 


3.42 


The  solution  will  lie  vrithin  an  exponentially  decreasing  envelope  vdiich 
has  a  time  constant  of  l/(  ) .  This  danped  oscillation  is  shown  in  Figure 
3.13.  From  Equation  3.91  and  Figure  3.14,  note  that  the  numerical  value  of 
danping  ratio  has  a  powerful  effect  on  system  response. 


FIGURE  3.13.  SECOSID  ORDER  DAMPED  OSCILLATIONS 


If  Equation  3.81  is  divided  M 


X 


X 


f(t) 

M 


(3.92) 


or,  rewriting  using  co^  and  Z,  defined  by  Equations  3.85  and  3.86 

V+2Cw„x  +  &)^x  =  (3.93) 

Equation  3.93  is  a  form  of  Equation  3.81  that  is  useful  in  analyzing  the 
behavior  of  any  second  order  linear  system.  In  general,  the  magnitude  and 
sign  of  damping  ratio  determine  the  response  properties  of  the  system.  There 
are  five  distinct  cases  \diich  are  given  names  descriptive  of  the  response 
associated  with  each  case.  These  are: 
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1.  C  =  0, 

undanped  A 

2.  0  <  C  >  1, 

underdanped  ^ 

3.  C  =  1, 

critically  damped 

4.  C  >  1, 

overdanped 

5.  C  <  0, 

mstable. 

will  be  examined 

in  turn,  making  use  of  Equation  3.88,  repeated 

below 


®i,2  =  +  j«„  V 1  -  C"*  (3.88) 

3. 5. 3.1  Case  1;  t,  -  0,  Undamped.  For  this  condition,  the  roots  of  the 

characteristic  equation  are 

“l.2  =  ±  3% 


giving  a  transient  solution  of  the  form 


or 


X  =  c,  cos  0)  t  +  c,  sin  w  t 

t  1  n  2  n 


=  A  sin  (w^t  +  ♦) 


(3.94) 


(3.95) 


showing  the  system  to  have  the  transient  response  of  an  undanped  sinusoidal 
oscillation  with  frequency  oo^.  Hence,  the  designation  of  as  the  "undanped 
natural  frequency."  Figure  3.9  shows  an  undanped  system. 


3. 5. 3. 2  Case  2;  0  <  C<  1.0»  Underdamped.  For  this  case,  m  is  given  by 

Equation  3.88. 


(3.88) 
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Uie  transient  solution  has  the  form 

-Cto)  t  r - 

-  A  e  "  sin{ft)„  V  1  -  t  +  ♦)  (3.96) 

Ihis  solution  shows  that  the  system  oscillates  at  the  danped  frequency, 
and  is  bounded  by  an  exponentially  decreasing  envelope  with  time  constant 
l/(Cw„).  Figure  3.14  shows  the  effect  of  increasing  the  danping  ratio  from 


3. 5. 3. 3  Case  3;  C  =  l.Oy  Critically  Dan^d.  For  this  condition,  the 
of  the  characteristic  equation  are 

”®i,2  “  “n  (3.97) 

\diich  gives  a  transient  solution  of  the  form 

-CO  t  -CO  t 

Xt  =  c^e  "  +  Cjte  "  (3.98) 


This  is  called  the  critically  danped  case  and  generally  will  not  overshoot. 
It  should  be  noted,  however,  that  large  initial  values  of  x  can  cause  one 
overshoot.  Figure  3.14  shows  a  response  vdien  C  =  1.0. 

3. 5. 3. 4  Case  4;  C  >  1.0/  Overdanped.  In  this  case,  the  characteristic  roots 
are 


% 


2 


(3.99) 


which  shows  that  both  roots  are  real  and  negative.  The  system  will  have  a 
transient  vdiich  has  an  exponential  decay  without  sinusoidal  motion.  The 
transient  response  is  given  by 


This  response  can  also  be  written  as 

x^  =  c^e"^^"^!  +  (3.101) 

\diere  and  Xj  are  time  constants  for  each  exponential  term. 

This  solution  is  the  sum  of  two  decreasing  exponentials,  one  with  time 
constant  x^  and  the  other  with  time  constant  .  The  smaller  the  value  of  x, 


the  quicker  the  transient  decays.  Usually  the  larger  the  value  of  the 
larger  is  conpared  to  Tj  .  Figure  3.5  shows  an  overdaitped  system. 

3. 5. 3. 5  Case  5;  -  1.0  <  C  <  0/  Unstable.  For  the  first  Case  5  exanple,  the 

roots  of  the  characteristic  equation  are 


VT^  (3.102) 

These  roots  are  the  same  as  for  the  underdanped  case,  except  that  the 
exponential  term  in  the  transient  solution  shows  an  exponential  increase  with 
time. 


=  e  (c^  cos  t  +  Cj  sin  1  -  t)  (3.103) 

Whenever  a  term  appearing  in  the  transient  solution  grows  with  time  (and 
especially  an  exponential  growth),  the  system  is  generally  unstable.  This 
means  that  whenever  the  system  is  disturbed  from  equilibrium  the  disturbance 
will  increase  with  time.  Figure  3.11  shows  an  unstable  system 

Case  5;  C  =  -  1.0,  Unstable .  For  this  second  Case  5  exanple,  the 

roots  of  the  characteristic  equation  are 


”V,2  =+“n 

and 


CO  t 


\  =  e  "  (c^  +  Cjt) 

This  case  diverges  much  the  same  way  as  shown  in  Figure  3.7. 


(3.104) 


(3.105) 


(3.99) 
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The  response  can  be  written  as  the  sum  of  two  exponential  terms 

m  t  m  t 

x^.  =  ^  +  tc^e  ^ 

where  the  values  of  m  can  be  determined  from  Equation  3.99. 

Case  5:  C  <  -  1.0>  Unstable.  This  third  Case  5  exanple  is  similar  to 
Case  4,  except  that  the  system  diverges  as  shown  in  Figure  3.7. 

®1,2  =  -  1  (3.99) 

The  response  can  be  written  as  the  sum  of  two  exponential  terms 

m  t  m  t 
Xt  =  c,e  "  +  c^e  ' 

vAiere  the  values  of  m  can  be  determined  from  Equation  3.99. 

Five  exanples  will  illustrate  some  of  these  system  response  cases. 

EXAMPLE 

Given  the  homogeneous  equation, 

•  • 

X  +  4x  =  0 


from  Equation  3.93, 


and 


C  =  0 


cci  =  2.0 

n 
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The  system  is  undan^d  with  a  solution 


=  A  sin(2t  +  «f>) 

\diere  A  and  are  constants  of  integration  which  could  be  determined  by 
substituting  boundary  conditions  into  the  total  general  solution. 

EXAMPLE 


Given  the  homogeneous  equation 

•  •  • 

X  +  X  +  X  =  0 

from  Equation  3.93, 

“n  =  1-0 

and 

C  =  0.5 

Also  from  Equation  3.87,  the  definition  of  damped  frequency 

%  =  (0^  V  1  -  =  0.87 

The  system  is  underdartped  with  a  solution 

x^  =  Ae  sin(0.87  t  + 

EXAMPLE 

Given  the  homogeneous  equation 

•  * 

I  +  X  +  X  =  0 

Multiply  by  four  to  get  the  equation  in  the  form  of  Equation  3.93. 
Then 

X  +  4x  +  4x  =  0 
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and 


w  =  2.0 

n 


C  =  1.0 

The  system  is  critically  danped  and  has  a  solution  given  by 


+  Cjte 


EXAMPLE 

Given  the  hcanogeneous  equation 


y  +  8x  +  4x  =  0 


from  Equation  3.93, 


and 


(0  ==  2.0 

n 


C  =  2.0 


The  system  is  overdanped  and  has  a  solution 


Xt  Cj^e  ^2® 


EXAMPLE 

Given  the  homogeneous  equation 


X  -  2x  +  4x  =  0 


from  Equation  3.93, 


0)  =  2.0 

n 


and 

C  =  -0.5 

From  Equation  3.87,  the  definition  of  damped  frequency 

=  1-7 

The  solution  is  unstable  (negative  damping)  and  has  the  form 

Xj.  =  Ae^  sin  (1.7  t  +  <|>) 
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In  summary,  the  best  danping  ratio  for  a  system  is  determined  by  the 
intended  use  of  the  system.  If  a  fast  response  is  desired  and  the  size  and 
nrunber  of  overshoots  is  inconsequential,  then  a  small  value  of  damping  ratio 
would  be  desired.  If  it  is  essential  that  the  system  not  overshoot  and 
response  time  is  not  too  critical,  a  critically  danped  (or  even  an  overdanped) 
system  could  be  used.  The  value  of  damping  ratio  of  0.7  is  often  referred  to 
as  an  optimum  damping  ratio  since  it  gives  a  small  overshoot  and  a  relatively 
c[uick  response.  The  optimum  damping  ratio  will  change  as  the  requirements  of 
the  physical  system  change. 

3.6  ANALOGOUS  SECCM)  ORDER  LINEAR  SYSTEMS 

3.6.1  Mechanical  System 

The  second  order  equation  \diich  has  been  examdned  in  detail  represents 
the  mass-spring-damper  system  of  Figure  3.12  and  has  a  differential  equation 
vdiich  was  given  by 


Ite  +  Dx  +  Kx  =  f(t) 


Using  the  definitions 


C 


- (3.85)  and  w 

2  \/m 


Equation  3.81  was  rewritten  as 

*  2 

X  +  2Cw  X  +  w  X 
n  n 


f(t) 

M 


(3.81) 


(3.86) 


(3.93) 


3.6.2  Electrical  System 

Hie  second  order  equation  can  also  be  applied  to  the  series  LRC  circuit 
shown  in  Figure  3.15. 


FIGURE  3.15.  SERIES  ELECTRICAL  CIRCUIT 


vdiere 

L  =  inductance 
R  =  resistance 
C  =  capacitance 
q  =  charge 
i  =  current 

Assume  q(0)  =  q(0)  =  0,  then  Kirchhoff's  voltage  law  gives 


Since 


‘  -  Pt 

(3.106) 


(3.107) 

C  s  - (3.108) 

2  VL/C 

and 

=  I  (3.109) 

Using  these  parameters.  Equation  3.106  can  be  written 

q  +  2Cw^q  +  q  =  (3.110) 


E(t)  =  Iq  +  Rq  +  3 

The  following  parameters  can  now  be  defined 


3.6.3  Se  rvomechani sms 

For  linear  control  systems  work  in  Chapter  13,  the  applicable  second 
order  equation  is 


(3.111) 

v^ere 

I  =  inertia 
f  =  friction 
/j  =  gain 
=  input 
©Q  =  output 
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Rearranging  Equation  3.111 


+  I  ®o  +  I  ®o  =  I  ®i  (3.112) 

or 

a  +  2C&)  a  +  0)  ^e.  =  w  ^a 

0  nO  nO  ni 

vrtiere  the  following  parameters  are  defined 


2  ^/7T 

Thus,  in  general,  any  second  order  differential  equation  can  be  written  in  the 
form 


(3.113) 

(3.114) 

(3.115) 


X  +  2?;a)^x  +  =  fj(t)  (3.116) 

vAiere  each  term  has  the  same  qualitative  significance,  but  different  physical 
significance.  " 

3.7  lAPIACE  TRANSFORMS 

A  technique  has  been  presented  for  solving  linear  differential  equations 
with  constant  coefficients,  with  and  without  inputs  or  forcing  functions.  The 
method  has  limitations.  It  is  suited  for  differential  equations  with  inputs 
of  only  certain  forms.  Further,  solution  procedures  require  looking  for 
special  cases  vdiich  require  careful  handling.  However,  these  procedures  have 
the  remarkable  property  of  changing  or  "transforming"  a  problem  of  integration 
into  a  problem  in  algebra,  that  is,  solving  a  quadratic  equation  in  the  case 
of  linear  second  order  differential  equations.  This  is  acconplished  by  making 
an  assunption  involving  the  number  e. 
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Given  the  second  order  homogeneous  equation 

••  • 

ax  +  bx  +  cx  =  0  (3.117) 

The  following  solution  is  assumed 

=  e"’^  (3.118) 

Substituting  into  Equation  3.117  gives 

am^e"^  +  bme"^  +  ce"^  =  0  (3.119) 

and,  factoring  the  exponential  term 

e”*(am^  +  Im  +  c)  =  0  (3.120) 

leading  to  the  assertion  that  Equation  3.118  will  produce  a  solution  to 
Equation  3.117  if  m  is  a  root  of  the  characteristic  equation 

am^  +  bm  +  c  =  0  (3.121) 

Introducing  operator  notation,  p  =  d/dt,  the  characteristic  equation  can  be 
written  by  inspection. 

ap^  +  bp  +  c  =0  (3.122) 

Equation  3.122  can  then  be  solved  for  p  to  give  a  solution  of  the  form 

Pit  p  t 

x^  =  Cie  "  +  c^e  (3.123) 

Of  course,  the  great  shortcoming  of  this  method  is  that  it  does  not  provide  a 
solution  to  an  equation  of  the  form 

ax  +  bx  +  cx  =  f(t)  (3.124,) 
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It  works  only  for  the  homogeneous  equation.  Still,  a  solution  to  the  equation 
can  be  found  by  obtaining  a  particular  solution  and  adding  it  to  the  transient 
solution  of  the  homogeneous  equation.  The  technique  used  to  obtain  the 
particular  solution,  the  method  of  undetermined  coefficients,  also  provides  a 
solution  by  algebraic  manipulation. 

However,  there  is  a  technique  v^ich  exchanges  (transforms)  the  \^ole 
differential  equation,  including  the  input  and  initial  conditions  into  an 
algebra  problem.  Fortunately,  the  method  applies  to  linear  first  and  second 
order  equations  with  constant  coefficients. 

In  Equation  3.124,  x  is  a  function  of  t.  For  enphasis.  Equation  3.124 
can  be  rewritten 


ax(t)  +  bx(t)  +  cx(t)  *  f(t) 


(3.125) 


Multiplying  each  term  of  Equation  3.125  by  the  integrating  factor  e"^  gives 


^(t)e”'^  +  bx(t)e'*’^  +  cx(t)e"’^ 


=  f(t)e' 


n  t 


(3.126) 


It  is  now  possible  that  Equation  3.126  can  be  integrated  term  by  term  on  botl^^B 
sides  of  the  equation  to  produce  an  algebraic  expression  in  m.  The  algebraic 
expression  can  then  be  manipulated  to  eventually  obtain  the  solution  of 
Equation  3.125. 

The  new  integrating  factor  e”'"  should  be  distinguished  from  the  previous 
integrating  factor  used  in  developing  the  operator  techniques  for  solving  the 
homogeneous  equation.  In  order  to  acconplish  this,  m  will  be  replaced  by  -s. 

The  reason  for  the  minus  sign  will  be  apparent  later.  In  order  to  integrate 
the  terms  in  Equation  3.126,  limits  of  integration  are  required.  In  most 
physical  problems,  events  of  interest  take  place  siibsequent  to  a  given 
starting  time  vhich  is  called  t  =  0.  To  be  sure  to  include  the  duration  of 
all  significant  events,  the  conposite  of  effects  from  time  t  =  0  to  time  t  =  «> 
will  be  included.  Equation  3.126  now  becomes 


% 

ax(t)  e“®^  dt  + 

>  p 

b  x(t)  e~®^  dt  + 

W  h 

c  x(t)  e 


-St  dt 
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dt 


(3.127) 


■» 

=  f(t) 

0 

Ecjuation  3.127  is  called  the  Laplace  transform  of  Ec[uation  3.125.  The 
problem  now  is  to  integrate  the  terms  in  the  equation. 

3.7.1  Finding  the  Laplace  Transform  of  a  Differential  Equation 

The  integrals  of  the  terms  of  Equation  3.127  must  now  be  found.  The 
Laplace  transform  is  defined  as 


€0 

x(t)  e"®’'  dt  =  L{x(t)}  =  X(s)  (3.128) 

V 

0 


vdiere  the  letter  L  is  used  to  signify  a  Laplace  transform.  X(s)  must,  for  the 
present,  remain  an  unknown.  (m  was  carried  along  as  an  unknown  until  the 
characteristic  equation  evolved,  at  vdiich  time  m  was  solved  for  explicitly. ) 
Since  Equation  3.128  transforms  x(t)  into  a  function  of  the  variable,  s. 


then 


OP 


OP 


c  x(t)  e  dt 

0 


c 


x(t)  dt 


0 


cX(s) 


(3.129) 


and  X(s)  will  be  carried  along  mtil  such  time  that  it  can  be  solved  for. 
The  transform  for  the  second  term,  bx(t)  is  given  by 


00 


00 


J  b  x(t)  dt 
0 


b 


x(t)  e"^’^ 

0 


dt 


(3.130) 
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To  solve  Equation  3.130,  a  useful  fonnula  known  as  integration  ty  part 
is  used 


p 

udv  =  rtv 

- 

w 

vdu 


(3.131) 


implying  this  formula  to  Equation  3.130,  let 


u  =  e 


St 


and 


then 


and 


dv  =  x(t)  dt 


du 


=  -se"®^  dt 


V  =  x(t) 

Substituting  these  values  into  Equation  3.131  and  integrating  from  t  =  0  to 

t  =  «o 


00  00 


x(t)e"®^  dt 


=  x(t)e' 


St 


x(t)(-se"®''  )dt 


=  x(t)e 


-St 


00  00 

+  S 
0  0 


x(t)e"®’^  dt 
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00 


=  x(t)e' 


■St 


+  sX(s) 


Now 


(3.132) 


lim  x(t)e“®^  -  x(0) 
t  ->  <» 


(3.133) 


and  assume  that  the  term  e"®'^  "dominates"  the  term  x(t)  as  t  ®.  The  reason 
for  using  the  minus  sign  in  the  exponent  should  now  be  apparent.  Thus, 
lim  x(t)e  =  0,  and  Equation  3.131  laecomes 
t  ^  ® 


x(t)e"®’'  dt  =  0  -  x(0)  +  sX(s)  =  sX(s)  -  x(0)  (3.134) 

% 

0 

Equations  3.129  and  3.134  can  now  be  abbreviated  to  signify  Laplace 
transformations . 


L{x(t)} 

=  X(s) 

(3.135) 

L{cx(t)} 

=  cX(s) 

(3.136) 

L{x(t)} 

=  sX(s)  -  x(0) 

(3.137) 

L{bx(t)} 

=  bfsX(s)  -  x(0)'l 

(3.138) 

Equation  3.138  can  be  extended  to  higher  order  derivatives.  Such  an  extension 
gives 

L{ax(t)}  =  a  [s^X(s)  -  sx(0)  -  x(0)]  (3.139) 
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Returning  to  Equation  3.127,  note  that  the  Laplace  transforms  of  all  the  terma 
except  the  forcing  function  have  been  found.  To  solve  this  transform,  thl 
forcing  function  must  be  specified.  A  few  typical  forcing  factions  will  be 
considered  to  illustrate  the  technique  for  finding  Laplace  transforms. 


EXAMPLE 


Then 


f(t)  =  A  »  constant 


L{A}  =  .  Ae"®^  dt  =--  e'^'^i-sdt) 


A  -St 
—  e 


L{A}  =  f 


(3.140) 


EXAMPLE 


f(t)  =  t 


Then 


L{t}  =  te"  dt 


To  integrate  by  parts,  let 


Then 


u  =  t 


dv  =  e”  dt 


du  =  dt 


1  -St 

V  =  --e 
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Substituting  into  Equation  3.131 


or 


EXAMPLE 


Then 


or 


EXAMPLE 


Then 


J  te"®^  dt 
0 


L{t}  =  ^ 
s 

f(t)  =  e""^ 


(3.141) 


Me""} 


2t  -St 

e  e  dt 


2-s  )  t 


dt 


1 


s  -  2 


Me""} 


(3.142) 


f(t)  =  sin  at 


L{sin  at} 


00 

sin  at  e“®"  dt 


0 
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Integrate  parts,  letting 


u  =  sin  at 


dv  =  e~°’^  dt 


Then 


du  =  a  (cos  at)dt 


1  -St 

V  =  -i  e 


Substituting  into  Ekjuation  3.131 


(sin  t)  e' 


■St  ^  -(sin  at)(e 
s 


^  I 


(cos  at)  e”®’'  dt 


or 


(sin  at)  e  dt  =  0  +  f 

5 


(cos  at)  dt  ( 


The  expression  (cos  at)  can  also  be  integrated  by  parts,  letting 


and 


u  =  cos  at 


dv  =  e’®^  dt 


du  =  -a  (sin  at)dt 


V  =  i  e-' 


3.143) 
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Giving 


{nr.c.  at)  dt  - 

-  a 

s 

s 

0  c 

0 

(sin  at)  e  dt 


or 


(cos  at)  e  dt  = 


I  -  I  L{sin  at) 


(3.144) 


Substituting  Equation  3.144  into  Equation  3.143 


L{sin  at}  =  0  +  I  (  i  -  I  L{sin  at})  = 

o  S  S 


s^  ~  s^ 


L{sin  at} 


vdiich  "obviously"  yields 


L{sin  at}  =  — -  (3.145) 

s  +  a 

Also  note  that  Equation  3.143  may  be  written  as 

L{sin  at}  =  -  L{cos  at} 
s 

vdiich  yields 

L{cos  at}  =  ^  (3.146) 

s  +  a 

The  Laplace  transforms  of  more  coirplicated  functions  may  be  quite  tedious 
to  derive,  but  the  procedure  is  similar  to  that  above.  Fortunately,  it  is  not 
necessary  to  derive  Laplace  transforms  each  time  they  are  needed.  Extensive 
tables  of  transforms  exist  in  most  advcinced  roathematics  and  control  system 
textbooks.  All  of  the  transforms  needed  for  this  course  are  listed  in  Table 
3.2  Page  3.73. 

The  technique  of  using  Laplace  transforms  to  assist  in  the  solution  of  a 
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differential  equation  is  best  described  by  an  example. 


EXAMPLE 


Given  the  differential  equation 

x*  +  4x  +  4x  =  4e“  (3.147) 

with  initial  conditions  x(0)  =  1,  x(0)  =  -4.  Taking  the  Laplace  transform 

of  the  equation  gives 

s^X(s)  -  sx(0)  -  x(0)  +  4  IsX(s)  -  x(0)]  +  4x(s)  =  — ^ 


or 


[s^  4  4s  +  4]  X(s)  4  [-S  4  4  -  4]  = 


s  -  2 


Solving  for  X(s) 


X(s)  = 


s  -  2s  4  4 
(s  -  2)(s  4  2y 


(3.148) 


In  order  to  continue  with  the  solution,  it  is  necessary  to  discusj 
partial  fraction  expansions. 


3.7.2  Partial  Fractions 

The  method  of  partial  fractions  enables  the  separation  of  a  complicated 
rational  proper  fraction  into  a  sum  of  simpler  fractions.  If  the  fraction  is 
not  proper  (the  degree  of  the  mmnerator  less  than  the  degree  of  the 
denominator),  it  can  be  made  proper  by  dividing  the  fraction  and  considering 
the  remaining  expression.  Given  a  fraction  of  two  polynomials  in  the  variable 
s  as  shown  in  Equation  3.148  there  occur  several  cases; 

3. 7. 2.1  Case  1;  Distinct  Linear  Factors.  To  each  linear  factor  such  as 
(as  4  b),  occurring  once  in  the  denominator,  there  corresponds  a  single 
partial  fraction  of  the  form,  Pi/{as  4  b). 
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EXAMPLE 


?rs  -M(s  -H  2)  =  (3.149) 
vAiere  A,  B,  and  C  are  constants  to  be  determined. 

3. 7. 2. 2  Case  2;  Repeated  Linear  Factors.  To  each  linear  factor,  (as  +  b), 
occurring  n  times  in  the  denominator  there  corresponds  a  set  of  n  partial 
fractions. 


EXAMPLE 

S^  -  9s  17  ^  A  ^  B  ^  C 

(s  -  2)^(s  +  1)  (s  +  1)  (s  -  2)  (s  -  2)^ 


(3.150) 


v4iere  A,  B,  and  C  are  constants  to  be  determined. 

3. 7. 2. 3  Case  3;  Distinct  Quadratic  Factors.  To  each  irreducible  quadratic 
factor,  as  +  bs  +  c,  occurring  once  in  the  denominator,  there  corresponds  a 
single  partial  fraction  of  the  form,  (As  +  B)/(as^  +  bs  +  c). 


EXAMPLE 

3s^  +  5s  +  8  A  Bs  +  C 

- i -  =  -  +-7 - ^  (3.151) 

(s  +  2)(s  +1)  s  +  2  s  +  1 

v4iere  A,  B,  and  C  are  constants  to  be  determined. 

3 *7. 2. 4  Case  4;  Repeated  Quadratic  Factors.  To  each  irreducible  quadratic 
factor,  as  +  bs  +  c,  occurring  n  times  in  the  denominator,  there  corresponds 
a  set  of  n  partial  fractions. 


EXAMPLE 

10  s^  +  s  4-  36  ^  A  Bs  4  C  ^  Ds  4  E 

(s  -  4)(s^  4  4)^  s  -  4  s^  4  4  (s^  4  if 

v^iere  A,  B,  C,  D,  and  E  are  constants  to  be  determined. 
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The  "brute-force"  technique  for  finding  the  constants  will  be  illustrab 
by  solving  Equation  3.152.  Start  by  finding  the  common  denominator  on  tl 
right  side  of  Equation  3.152 


10  S^  S  +  36  ^  A(s^  +  4)^  +  (Bs  +  C)(s  -  4)(s^  +  4)  (Ds  +  E)(s  -  4) 

(s-4)(s'+4)'  (s-4)(s^+4)"  (3.153) 

Then  the  numerators  are  set  equal  to  each  other 

lOs^  +  s  +  36  =  A(s^  +  4)^  +  (Bs  +  C)(s^  +  4)(s  -  4)  +  (Ds  +  E)(s  -  4) 

(3.154) 

Since  Equation  3.154  must  hold  for  all  values  of  s,  enough  values  of  s  are 
substituted  into  Equation  3.154  to  find  the  five  constants. 

1.  Let  s  =  4,  then  Equation  3.154  becomes 


and 


(10)(16)  +4+36 


=  400A 


A  =  1/2 

2.  Let  s  =  2j,  then  Equation  3.154  becomes 

-40  +  2j  +  36  =  -4D  +  2jE  -  8jD  -  4E 
-  4  +  2j  =  -4(D  +  E)  +  2j  (E  -4D) 

Ihe  real  and  imaginary  parts  must  be  equal  to  their  counterparts  on  the 
opposite  side  of  the  equal  sign,  thus 


and 


(D  +  E)  =  1 

E  -  4D  =  1 
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or 

D  =  0 

and 

E  =  1 

3.  Now  let  s  =  0,  then  Equation  3.154  becomes 

36  =  16A  -  16  (C)  -  4E 

and  from  steps  1  and  2 

A  =  1/2,  E  =  1 

hence 

36  =  8  -  16C  -  4 

and 

C  =  -2 

4.  Let  s  =  1,  then  Equation  3.154  becomes 

47  =  25  (1/2)  +  (B  -  2)  (-15)  -3 

94  =  25  -  30B  +  60-6 
or 

B  =  -  1/2 

Now  Equation  3.155  may  be  written  by  siibstituting  the  values  of  A, 
and  E  into  Equation  3.152 


10  s^  +  s  +  36  _  1  -1/2  s  +  4  ^  1 

(s  -4)(s^  +4)"  s  -  4  s'  +  4  (s'  +  4)' 


B,  C,  D, 


(3.155) 
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Returning  now  to  the  example  Laplace  solution  of  the  differential  equation 

(3.147) 


X  +  4x  +  4x  =  4e^^ 


Ihe  Laplace  transformed  equation  was 


X(s)  = 


s  -  2s  4 
(s  -  2)(s  +  2)^ 


(3.148) 


which  can  now  be  expanded  by  partial  fractions 


s^  -  2s  +  4  ^  A  ^  B  ^  C 

(s  -  2)(s  +2)^  s  -  2  s  +  2  (s  +  2)^ 


(3.156) 


Taking  the  common  denominator,  and  setting  numerators  equal 

s^  -  2s  +  4  =  A(s  +  2)^  +  B(s  +  2)(s  -  2)  +  C(s  -  2)  (3.157) 

The  "brute-force"  technique  could  again  be  used  to  solve  for  the 
constants  A,  B,  and  C  by  substituting  different  values  of  s  into  Equation 
3.157.  An  alternate  method  exists  for  solving  for  the  constants.  Multiplying 
the  right  side  of  Equation  3.157  gives 

S^  -  2s  +  4  =  As^  +  4As  +  4A  +  Bs^  -  4B  +  Cs  -  2C  ^ 


S^  -  2s  +  4  =  (A  +  B)s^  +  (4A  +  C)S  +  (4A  -  4B  -  2C) 


Now  the  coefficients  of  like  powers  of  s  on  both  sides  of  the  equation  must  be 
equal  (that  is,  the  coefficient  of  s^  on  the  left  side  equals  the  coefficient 
of  s^  on  the  right  side,  etc.).  Equating  gives 

s^  ;  1  =  A  +  B 

s^  :  -2  =  4A  +  C 

S°  ;  4  =  4A  -  4B  -  2C 
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Solving  for  the  constants  gives 


A  =  1/4 
B  =  3/4 

C  =  -3 

Substituting  the  constants  into  Equation  3.156  results  in  the  expanded  right 
side 


X(s) 


1/4 


1 

s^ 


+  3/4 


1 

s  +  2 


-3 


1 

s  +  2 


2 


(3.158) 


Another  expansion  method  called  the  Heaviside  Expansion  Theorem  can  be  used  to 
solve  for  the  constants  in  the  numerator  of  distinct  linear  factors.  This 
method  of  expansion  is  used  extensively  in  Chapter  13,  Linear  Control  Theory. 
If  the  denominator  of  an  expansion  term  has  a  distinct  linear  factor,  (s  -  a), 
the  constant  for  that  factor  can  be  found  by  multiplying  X(s)  by  (s  -  a)  and 
evaluating  the  remainder  of  X(s)  at  s  =  a. 

Stated  mathematically  the  Heaviside  Expansion  Theorem  is 


A  =  (s  -  a)  X(s) 


s  =  a 


EXAMPLE 

X(s) 


7s  -  4 


s(s  -  l)(s  +2)  s  s-1  s+2 


A  =  sX(s) 


7s  -  4 


-4 


s  =  0 


(s  -  l)(s  +2)  ^  ^  ^  -  (-1)(2) 


=  2 
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B  =  (S  -  I)X(S) 


Is  -  A  _  1  -  A  _ 

S(S  +  2)  3,1"  TUU)  - 


C  »  (s  +  2)X(s) 

s  =  -2 


7s  -  4 
s(s-  1) 


-14  -  4  _  , 

“  (-2) (-3)  = 


As  another  exanple,  the  constant  A  in  the  first  term  on  the  right  side  of 
Equation  3.156  can  be  evaluated  using  the  Heaviside  Expansion  Theorem. 


s^  -  2s  +  4  ^  A  ^  B  ^  C 

(s  -  2)(s  +  2)^  s  -  2  s  +  2  (s  +  2)^ 


(3.156) 


A  =  (s  -  2)X(s)  ^  s^  -  2s  -h  4  ^4  ^1 

s  =  2  is  +  2f  s  =  2  ^  ^ 

>Aiich  is  the  same  result  obtained  earlier  by  equating  like  powers  of  s. 

3.7.3  Finding  the  Inverse  Laplace  Transform  ^ 

Now  that  methods  to  expand  the  right  side  of  X(s)  have  been  discussed  in 
detail,  all  that  remains  is  to  transform  the  expanded  terms  back  to  the  time 
domain.  This  is  easily  acconplished  using  any  suitable  transform  table. 

Returning  to  the  Laplace  transformed  and  expanded  equation  in  the  exanple 


X(s)  = 

Using  Table  3.2, 
transformed  to 


1/4 


s  -  2 


+  3/4 


s  +  2 


-  3 


(s  +  2) 


(3.158) 


it  can  be  easily  verified  that  Equation  3.158  can  be 


x(t)  =  1/4  e^’^  +  3/4  e~“  -  3t  e~“ 


(3.159) 
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In  svnnmary,  the  strength  of  the  Laplace  transform  is  that  it  converts 
linear  differential  equations  with  constant  coefficients  into  algebraic 
equations  in  the  s-domain.  All  that  remains  to  do  is  to  take  the  inverse 
transform  of  the  explicit  solutions  to  return  to  the  time  domain.  Although 
the  applications  here  at  the  School  will  consider  time  as  the  independent 
variable,  a  linear  differential  equation  with  any  independent  variable  may  be 
solved  by  Laplace  transforms. 

3.7.4  Laplace  Transform  Properties 

There  are  several  inportant  properties  of  the  Laplace  transform  \(diich 
should  be  included  in  this  discussion. 

In  the  general  case 


=  s^  x(s)  - 

'  dt“  ; 


[s*"'  x(0)  +  + 


,  d  ^  ^x(0)j 


dt' 


n-l 


(3.160) 


For  quiescent  systems 


(3.161) 


This  result  enables  transfer  functions  to  be  written  by  inspection. 
EXAMPLE 

Given  the  differential  equation 

*x*  +  4x  +  4x  =  40^*"  (3.162) 

with  quiescent  initial  conditions,  the  Laplace  transform  can  immediately  be 
written  by  inspection  as 

X(s)(s^  +  4s  +  4)  =  (3.163) 
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In  roost  cases,  reference  to  Table  3.2  will  probably  be  needed  to  transform 
right  side  forcing  function  (input). 

Another  significant  transform  is  that  of  an  indefinite  integral.  In  the 
general  case 


..  x(t)dt" 


) 


n-1 


+..(3.164) 


Equation  3.164  allows  the  transformation  of  integro-differential  equations 
such  as  those  arising  in  electrical  engineering. 

For  the  case  vdiere  all  integrals  of  f(t)  evaluated  at  zero  are  zero 
(quiescent  system)  the  transform  becomes 


jjll  -  ^ 


EXAMPLE 


Given  the  differential  equation 


« #  •  r  2 1 

X  +  4x  +  4x  +  xdt  =  4e 


(3.165) 


(3.166) 
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TABLE  3.2. 
LAPIACE  TRANSFORMS 


1.  aX(s) 


ax{t) 


2.  a[sX(s)  -x(0)] 


ax(t) 


3.  a[s‘^(s)  -  sx(0)  -  x(0)] 


ax(t) 


(which  can  be  extended  to 
ary  necessary  order) 


(n  =  1,  2, 


s  +  a 


(s  +  a)' 


(s  +  a) 


(n  =  1,  2,  ...) 


t^  e-^"^ 


(s  +  a)  (s  +  b)  ^  ^ 


b^  (e-^^  -  e 


(s  +  a)  (s  +  b) 


a  #  b 


^  (ae-^^  -  be-^tj 


(s  +  a)  (s  +  b)  (s  +  c) 


(b-c)e  -  (a-c)e  +  (a-b)e 
(a-b) (b-c) (a-c) 
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TABLE  3.2 

LAPLACE  TRANSFORMS  (continued) 
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with  quiescent  initial  conditions,  the  Laplace  transform  can  inrmediately  be 
written  by  inspection  as 


X(s)  (s^  +  4s  +  4)  +  ~  —1^5 

The  right  side  transform  is  the  same  as  Equation  3.163.  Factoring  results  in 

X(s)  (s^+4s  +  4+|)  =  (3.167) 

Multiplying  Equation  3.167  ty  s  gives 

X(s)  (s^  +  4s^  +  4s  +  1)  =  — ^ 


which  raises  the  order  of  the  left  side  and  acts  to  differentiate  the  right 
side  (input). 

The  usefulness  of  the  Laplace  transform  technique  will  be  demonstrated  by 
solving  several  example  problems. 


EXAMPLE 


Solve  the  given  equation  for  x(t), 

X  +  2x  =  1  (3.168) 

\dien  x(0)  =  1. 

By  Laplace  transform  of  Equation  3.168 

L{x}  =  sX(s)  -  x(0) 

L{2x}  =  2X(s) 

L(l)  =  i 
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Thus 


(s  +  2)  X{s)  =  ^  +  1 

5 


Y/e\  —  S  +  1  _  A  B 

"  s(s  +  2)  “  s  ^  s  +  2 


Solving, 


and 


A  =  1/2 


B  =  1  -  1/2  =  1/2 


Inverse  Laplace  transforming  gives 


x{t)  «  1/2  4  1/2  e' 


-2t 


(3.169) 


EXAMPLE 


Given  the  differential  equation 


X  +  2x  =  sin  t,  x(0)  =  5 


(3.170) 


solve  for  x(t). 


Taking  the  Laplace  transform  of  Equation  3.170 


and 


sX(s)  -  x(0)  +  2X(s)  = 


X(s)  = 


s^  +  1 


(s^  +  l)(s  +2)  s  +  2 


(3.171) 
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Expanding  the  first  term  on  the  right  side  of  the  equation  gives 


_ 1 _  ^  As  B  ^  C 

(s^  +  l)(s  +2)  s^  +  1  s  +  2 


Taking  the  common  denominator  and  equating  nmerators  gives 

1  =  (As  +  B)(s  +  2)  +  C(s^  +  1) 
Substituting  values  of  s  leads  to 


A  =  -  1/5 
B  =  2/5 
C  =  1/5 


and  sijbstituting  back  into  Equation  3.171  gives 

X(S)  .  .  _1ZS_  ^  5 


s  +1  s^+1  s  +  2  s  +  2 


Inverse  Laplace  transforming  gives  the  solution 

X(t)  =  -1/5  cos  t  +  2/5  sin  t  +  5  1/5  e"“ 


EXAMPLE 

Given  the  differential  equation 

X  +  5x  +  6x  =  3e"^’^,  x(0)  =  x(0)  =  1 

solve  for  x(t). 

Taking  the  Laplace  transform  of  Equation  3.174 


s^X(s)  -  sx(0)  -  x(0)  +  5sX(s)  -  5x(0)  +  6X(s) 


3 

s  +  3 


(3.172) 


(3.173) 


(3.174) 


(3.175) 
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or 


X(s) 


+  9s  +  21 

(s  +  3)  (s^  +  5s  +  6) 


Factoring  the  denominator. 


X(s) 


_s^  +  9s  +  21__ 
(s  +  3){s  +  2)(s  +  3) 


(3.176) 


(3.177) 


X(s)  =  -  (3.178) 

(s  +  3)^s  +  2) 

X(s)  =— - ®_+— (3.179) 

s  +  3  (s  +  3)  s  +  2 

Finding  the  common  denominator  of  Ekjuation  3.179,  and  setting  the  resiiltant 
numerator  equal  to  the  numerator  of  Equation  3.178. 

s^  +  9s  +  21  =  A(s  +  3)(s  +  2)  +  B(s  +  2)  +  C(s  +  3)^ 

v^ich  can  be  solved  easily  for 

A  =  -6 

B  =  -3 

C  =  7 


s  +  3  (s  +  3)^  s  +  2 
vhich  can  be  inverse  Laplace  transformed  to 


Now  X(s)  is  given  by 

X(s)  = 
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x(t)  =  -  3te"^'^  +  7e"^^ 


(3.180) 


EXAMPLE 

Given  the  differential  equation 

*x*  +  2x  +  lOx  =  3t  +  6/10  (3.181) 

x(0)  =  3 

x(0)  =  -27/10 


solve  for  x(t). 


Laplace  transforming  Equation  3.181  and  solving  for  X(s)  gives 


X(s)  = 


3s  +  3.3s  +  0.6s  +  3 
s^  (s^  +  2s  +  10) 


^  +  CS  4  D 

S  S^  S^  4  2S  4  10 


(3.182) 


vdiere 


Thus, 


A 

B 

C 

D 


0 

0.3 

3 

3 


X(s)  =  -^  4 - 3sjtJ - 

S  S  4  2s  4  10 


(3.183) 


To  make  the  inverse  Laplace  transform  easier,  Equation  3.183  is  rewritten  as 

(s  4  1) 


X(s)  =  ^  +  3 


(S  4  1)^  4  3^ 


(3.184) 
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viliich  is  readily  inverse  transformable  to 


x(t)  =  0.3t  +  Se'^^cos  3t 

3.8  TRANSFER  FUNCTIOgS 


(3.185) 


A  transfer  function  is  defined  in  Chapter  13,  Linear  Control  Theory  as, 
"The  ratio  of  the  output  to  the  input  expressed  in  operator  or  Laplace 
notation  with  zero  initial  conditions."  The  term  "transfer  function"  can  be 
thought  of  as  vdiat  is  done  to  the  input  to  produce  the  output.  A  transfer 
function  is  essentially  a  mathematical  model  of  a  system  and  embodies  all  the 
physical  characteristics  of  the  system.  A  linear  system  can  be  conpletely 
described  by  its  transfer  fxjuiction.  Consider  the  following  quiescent  system. 

ax  +  bx  +  cx  =  f(t)  (3.186) 


x(0)  =  x(0)  =  0 


Talcing  the  Laplace  transform  of  Equation  3.186  results  in 

as^X(s)  +  bsx(s)  +  cx(s)  =  F(s) 

factoring  gives 

X(s)(as^  +  bs  +  c)  =  F(s) 


or 


X(s)  ^  1 

F(s)  as^  +  bs  +  c 


(3.187) 

(3.188) 


(3.189) 


Since  Equation  3.186  represents  a  system  vdiose  input  is  f(t)  and  whose  output 
is  x(t)  the  following  transforms  can  be  defined 

X(s)  s  output  transform 

F(s)  s  input  transform 
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The  transfer  function  can  then  be  given  the  symbol  TF  and  defined  as 


TF 


X(s) 

F(s) 


In  the  exanple  represented  by  Equation  3.189 


TF 


1 

as^  +  bs  +  c 


(3.190) 


(3.191) 


Note  that  the  denominator  of  the  transfer  function  is  algebraically  the  same 
as  the  characteristic  equation  appearing  in  the  Equation  3.186.  The 
characteristic  equation  conpletely  defines  the  transient  solution,  and  the 
total  solution  is  only  altered  by  the  effect  of  the  particular  solution  due  to 
the  input  (or  forcing  fxjnction).  Thus,  from  a  physical  standpoint,  the 
transfer  function  conpletely  characterizes  a  linear  system. 

Ihe  transfer  function  has  several  properties  that  will  be  used  in  control 
system  analysis.  Suppose  that  two  systems  are  characterized  by  the 
differential  equations 


ax*  +  bx  +  cx  =  f(t)  (3.192) 

and 

dy  +  ey  +  gy  =  x(t)  (3.193) 

From  the  equations  it  can  be  seen  that  the  first  system  has  an  input  f(t),  and 
an  output  x(t).  Ihe  second  system  has  an  input  x(t)  and  an  output  y(t).  Note 
that  the  input  to  the  second  system  is  the  output  of  the  first  system.  Taking 
the  Laplace  transform  of  these  two  equations  gives 

(as^  +  bs  +  c)  X(s)  =  F(s)  (3.194) 

and 

(ds^  +  es  +  g)  Y(s)  =  X(s)  (3.195) 
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Now,  both  of  these  systems  can  be  represented  schematically  by  the  block 
diagrams  shown  in  Figure  3.16. 


SYSTEM  2 


If  it  is  desired  to  find  the  output  y(t)  of  System  2  due  to  the  input 
f(t)  of  System  1,  it  is  not  necessary  to  find  x(t)  since  the  two  systems  can 
be  linked  using  transfer  functions  as  shown  in  Figure  3.17. 


FIGURE  3.17.  COMBINING  TRANSFER  FUNCTIC»FS 


The  solution  y(t)  is  then  given  by  the  inverse  transform  of  Y(s),  or 


Y(S)  = 

[TFj] 

F(S) 

(3.198) 

Y(S)  = 

ITFJ 

[TFj]  F(S) 

(3.199) 

This  method  of  solution  can  be  logically  extended  to  include  any  desired 
nijmber  of  systems. 

3.9  SIMULTftNBJUS  LINEAR  DIFFERENTIAL  EQUATIONS 

In  many  physical  problems  the  mathematical  description  of  the  system  can 
most  conveniently  be  written  as  simultaneous  differential  equations  with 
constant  coefficients.  The  basic  procedure  for  solving  a  system  of  n  ordinary 
differential  equations  in  n  dependent  variables  consists  in  obtaining  a  set  of 
equations  from  \diich  all  but  one  of  the  dependent  variables,  say  x,  can  be 
eliminated.  The  equation  resulting  from  the  elimination  is  then  solved  for 
the  variable  x.  Each  of  the  other  dependent  variables  is  then  obtained  in  a 
similar  manner. 

A  very  effective  means  of  handling  simultaneous  linear  differential 
equations  is  to  take  the  Laplace  transform  of  the  set  of  equations  and  reduce 
the  problem  to  a  set  of  algebraic  equations  that  can  be  solved  explicitly  for 
the  dependent  variable  in  s.  This  method  is  demonstrated  below. 

Given  the  set  of  equations 

j2  ,2 

3  — *  +  X  +  — y  +  3y  =  f(t)  (3.200) 

dt^  dt^ 

2^  +  x+^^2y  =  g(t)  (3.201) 

dV  dt^ 

v^ere  x(0)  =  x(0)  =  y(0)  =  y(0)  =  0,  find  x(t)  and  y(t).  Taking  the 
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Laplace  transform  of  this  system  yields 

(3s^  +  1)  X(s)  +  (s^  +  3)  Y(s)  =  F(s)  (3.202) 

(2s^  +  1)  X{s)  +  (s^  +  2)  Y(s)  =  G(s)  (3.203) 

Cramer's  nile  vdll  now  be  used  to  solve  this  set  of  equations.  Cramer's  rule 
can  be  stated  in  its  sinplest  form  as,  given  the  equations 


Pj(s)  X(s)  +  P2(s)  Y(s)  =  Fj(s) 
Qj(s)  X(s)  +  Q2(s)  Y(s)  =  F2(s) 


then, 


X(s)  = 


for  unknown  X(s),  and 


Y(s)  = 


for  the  xmknown  Y(s). 


(3.204) 

(3.205) 


(3.206) 


(3.207) 
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\diich  yields  a  solution 


x(t)  =  -2t  +  3  -  7/4e"^  -  l/4e^  +  1/2  sin  t  -  cos  t  (3.213) 


A  similar  approach  will  obtain  the  solution  for  y(t). 

In  the  case  of  three  simultaneous  differential  equations,  the  application 
of  Laplace  transforms  and  use  of  Cramer's  rule  will  yield  the  solution. 


Pj(s)X(s)  +  P2(s)Y(s)  +  P3(s)Z(s)  =  F^(s) 

Qj  (s)X(s)  +  Q2(s)Y(s)  +  Q3(s)z(s)  =  FjCs) 

R3(s)X(s)  +  R2(s)Y(s)  +  R3(s)Z{s)  =  F3(s) 

\diere 

Q3 

R3 

x(s)  =  - 

Px  ^2  P3 

Qx  Q2  Q3 

Rx  ^3 

Y(s)  and  Z(s)  will  have  similar  forms. 

3.10  ROOT  PLOTS 


(3.214) 

(3.215) 

(3.216) 


(3.217) 


Scrnie  insight  into  the  response  of  a  second  order  system  can  be  gained  by 
examining  the  roots  of  the  differential  equation  describing  the  system  on  a 
root  plot.  A  root  plot  is  a  plot  of  the  roots  of  the  characteristic  equation 
in  the  coitplex  plane.  Root  plots  are  used  in  Chapter  8,  Dynamics,  to  describe 
aircraft  longitudinal  and  lateral  directional  modes  of  motion.  These  plots 
are  also  used  extensively  in  Chapter  13,  Linear  Control  Theory,  for  linear 
control  system  analysis. 
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It  was  shown  earlier  that  a  second  order  linear  system  can  be  put  into 
the  following  form: 


X  +  2Cw  X  +  w^x  =  -4—  (3.93) 

^  n  n  M 

\rtiose  roots  can  be  written  as 

Pi. 2  =  -  ±  :j«n  (3.88) 

or 

Pi, 2  =  -  ±  (3.218) 

Figure  3.18  is  a  plot  of  the  two  roots  of  Equations  3.88  and  3.218  in  the 
conplex  plane. 


Pi. 


=  -  +  jw„  \[l^ 


(3.88) 


real 


2L, 


imaginary 


>rtiere  the  first  term  is  plotted  on  the  real  axis  and  the  second  term  plotted 
on  the  imaginary  (j)  axis. 


3.87 


3.88 


The  five  distinct  daitping  cases  previously  discussed  can  be  examined  on 
root  plots  through  the  use  of  Equation  3.88. 


1.  C  =  0,  undanped  (Figure  3.19) 

Pi. 2  =  -  ^  "n  ± 

Pi, 2  =  0  +  H 


Im 

Pi  •  An 


Re 


P2  •  -An 

I 


NEUTRALLY  STABLE 


(3.88) 


FIGURE  3.19.  NEUTRALLY  STABLE  UNDAMPED  RESPONSE 


2.  0  <  C  <  1.0,  Underdamped  (Figure  3.20) 

Pi, 2  =  -  ^  (3.88) 

Pi. 2  =  (-)  ±  j  (+) 
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X 


FIGURE  3.20.  STABLE  UNDERDAMPED  RESPONSE 


3.  C  =  1.0  Critically  danped  (Figure  3.21) 

Pi, 2  =  -  ±  j“„Vl  -  C" 

Pl  *5  =  -  CO 

*^1,2  n 


FIGURE  3.21.  STABLE  CRITICALLY  DAMPED  RESPONSE 
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-1.0  (Figure  3.23) 


X^  =  C,ennt  +  C2teAn‘  \ 

UNSTABLE 

FIGURE  3.23.  UNSTABLE  RESPONSE 
1.0  <  C  <  0  (Figure  3.24) 

,2  =  -  ±  3%^^  -  C' 

2  “  i  j 

,2  “  (+)  ±  j  (+) 


UNSTABLE 


FIGURE  3.24  UNSTABLE  RESPONSE 


C  <  -  1.0  Both  roots  positive  (Figure  3.25) 


Pi. 2  =  -  -  c"  (3.88) 

Pi. 2  =  (+) 


UNSTABLE 

FIGURE  3.25.  UNSTABLE  RESPONSE 


In  summary,  a  second  order  system  with  both  roots  located  to  the  left  of 
the  imaginary  axis  is  stable.  If  both  roots  are  on  the  imaginary  axis  the 
system  is  neutrally  stable,  and  if  one  or  more  roots  are  located  to  the  right 
of  the  imaginary  eixis  the  system  is  unstable.  These  conditions  are  shown  in 
Figure  3.26. 
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STABLE 


FIGURE  3.26.  ROOT  LOCUS  STABILITY 

Root  plots  can  be  used  for  analysis  of  the  aircraft  modes  of  motion.  For 
exanple,  the  longitudinal  static  statiblity  of  an  aircraft  is  greatly 
influenced  by  center  of  gravity  (eg)  position.  Figure  3.27  shows  how  the 
roots  of  the  characteristic  equation  describing  one  of  the  longitudinal  motion 
modes  change  position  as  the  eg  is  moved  aft.  This  plot  is  called  a  root 
locus  plot. 
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B 


f 

•  15%  MAC 


eg  LOCATION 

IMAGINARY 


•  20 


•  25 


REAL 


•  15%  MAC 


FIGURE  3.27  EFFECT  OF  CG  SHIFT  ON  LONGITUDINAL  STATIC 
STABILITY  OF  A  TYPICAL  AIRCRAFT 


Note  that  as  the  eg  is  moved  aft  of  its  initial  location  at -.15%  MAC,  danping 
of  this  mode  of  motion  (short  period)  increases  while  the  frequency  decreases. 
Zero  frequency  is  reached  between  a  eg  location  of  28%  and  30%  MAC.  The  root 
locus  then  splits  into  a  pair  of  real  roots,  branches  AB  and  AC  of  the  locus. 
These  branches  represent  danped  aperiodic  (nonoscillatory)  motion.  The  short 
period  mode  of  motion  goes  unstable  at  a  eg  location  of  35%  MAC.  The  location 
of  the  eg  vhere  this  instability  occurs  (35%  MAC  in  this  exanple)  is  known  as 
the  maneuver  point  and  it  is  discussed  in  detail  in  Chapter  6,  Maneuvering 
Flight. 
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PROBLEMS 


Solve  for  y. 


3.1.  .  X*  +  4x  +  sin  6x 


3.2.  — ^  »  e"*  +  sin  cox 
dx^ 

3.3  X* 

dx^ 

y  3x^  “  ® 

3.5.  (x  -1)^  ydx  +  x^  (y  -  l)dy  =  0 

Just  find  a  solution.  Solving  for  y  is  tough. 

Test  for  exactness  and  solve  if  exact. 

3.6.  (y^  -  X)  dx  +  (x^  -  y)  dy  =  0 

3.7.  (2x^  +  3y)  dx  +  (3x  +  y  -  1)  dy  =  0 

3.8.  (2xy"  e'^  +  2xy'  +  y)  dx  +  (x'y'e^  -  x'/  -  3x)  dy 

3.9.  Multiply  Problem  3.8  by  1/y  and  solve  for  y. 
y  ^  0. 


0 


Solve  for  y^ 


Note  this  assumes  that 


3.11.  Y'”  -  5y"  -  24y'  *  0 

3.12.  y"  +  12y'  +  36y  =  0 

3.13.  y"  +  4y'  +  13y  =  0 

Solve  for  y^  and  y^  in  Problems  3.14  -  3.17,  then  solve  for  the  general 


solution. 

3.14. 

•  •  •  5  ^ 

y  +  5y  +  6y  =  3e 

y(0)  = 

1,  y(0) 

=  6 

3.15. 

•  •  • 

y  +  4y  +  4y  =  cos  t 

y{0)  = 

S  '  y(0) 

104 

3.16. 

2x  +  4x  +  20x  =  6t  +  ^ 

,  X(0)  = 

3  ,  x(0) 

27 

■  "TO 

3.17. 

3x  +  2x  =  -4e~^^  , 

x(3)  = 

-0.14 

3.18. 

Find  co^ ,  co^ ,  C,  and 

describe  system 

dairping 

(i.e.,  underdamped 

overdanped,  etc.)  vdiere  applicable. 

••  •  _3t 

y  +  5y  +  6y  =  3e 

#*  • 

3.19.  y  +  4y  +  4y  =  cos  t 

3.20.  2x  +  4x  +  20x  =  6t  +  -g 

3.21.  3x  +  2x  =  -  4e'^*^ 
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In  problems  3.22  -  3.24  find  x(s),  do  not  find  the  inverse  transform. 


3.22. 

•  •  • 

3x  +  X  +  6x  = 

sin  6t, 

x{0) 

=  x(0)  = 

3.23. 

••  ■ 

X  -  2x  +  5x  = 

e”^  sin  3t, 

x(0) 

=  -1,  x(0) 

3.24. 

••  » 

4x  +  3x  -  X  = 

t^  -  t  sin  2t, 

x(0) 

=  3,  x(0) 

In  Problems  3.25  -  3.27,  expand  X(s)  by  partial  fractions  and  find  the 
inverse  transforms. 


3.25.  X(s)  = 


5s  +  29s  +  36 _ 

(s  +  2)  (s^  +  4s  +  3) 


3.26.  X(s)  = 


3.27.  X(s)  = 


2s  +  6s  +5 
(s^  +  3s  +  2)  (s  +  1) 

2s^  +  7s^  +  27s^  +  51s  +  27 
(s^  +  9s)  (s^  +  3s  +  3) 


Solve  the  following  problems  by  Laplace  transform  techniques. 


3.28.  X  +  2x  =  sin  t, 

•  •  • 

3.29.  X  +  5x  +  6x  =  3e 


•3t 


x(0)  =  5 

x(0)  =  x(0)  =  1 


Solve  using  Laplace  Transforms 


3.30.  X  +  3x  -  y  =  1 
X  +  8x  +  y  =  2 


x(0)  =  y(0)  =  0 
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3.31.  Read  the  question  and  circle  the  correct  answer,  true  (T)  or  false  (F): 


T 

T 


T 

T 

T 

T 

T 


T 

T 

T 

T 

T 


F 

F 


F 

F 

F 

F 

F 


F 

F 

F 

F 

F 


The  particular  solution  to  a  second  order  differential  equation 
contains  two  arbitrary  constants  that  are  solved  for  using 
initial  conditions  and  the  transient  solution  too. 

Solutions  to  linear  differential  equations  are  generally 
nonlinear  functions. 

Differential  equation  solutions  are  free  of  derivatives. 

Direct  integration  will  give  solutions  to  some  differential 
equations  without  the  necessity  of  arbitrary  constants. 

In  general,  the  number  of  arbitrary  constants  in  the  solution 
of  a  differential  equation  is  equal  to  the  order  of  the 
differential  equation. 

There  is  no  known  way  to  determine  if  a  differential  equation 
is  exact. 

The  solution  to  a  first  order  linear  differential  equation  with 
constant  coefficients  is  always  of  exponential  form. 

The  Laplace  variable  s  can  be  real,  imaginary,  or  conplex. 

Inverse  Laplace  transforms  are  used  to  return  from  the  s  to  the 
time  domain. 

Quiescent  systems  have  zero  initial  conditions. 

First  order  equation  roots  cannot  be  plotted  on  root  plots. 

A  transfer  function  can  be  defined  as  input  transform  divided 
by  output  transform. 

The  characteristic  equation  conpletely  describes  the  transient 
solution. 

The  method  of  mdetermined  coefficients  is  used  to  solve  for 
the  particular  solution. 

* 

X  +  4x  +  13x  =  3,  is  a  second  degree  equation, 

•#  • 

X  +  4x  +  13x  =  3,  is  a  second  order  equation. 

0 

4x  +  13x  =  3,  is  a  first  order  equation. 

It  is  inpossible  to  have  a  linear,  second  degree  equation. 

13x  =3,  is  a  linear  equation. 
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T  F  13x  *  3,  is  a  differential  equation.  ^ 

T  F  Daiqping  ratio  and  natural  frequency  have  no  physical 

significance. 

T  F  The  time  constant  and  time  to  half  anplitude  for  a  first  order 
system  are  equal. 

T  F  The  Laplace  transform  converts  a  differential  equation  from  the 
time  domain  to  the  s  domain. 


T 

T 

T 

T 

T 

T 

T 

T 

T 

T 

T 

T 

T 

T 

T 


F  The  transient  response  is  dependent  on  the  input. 

F  Laplace  transforms  are  easy  to  derive. 

F  In  general,  it  is  easier  to  check  a  candidate  solution  to  see 
if  it  is  a  solution  than  to  determine  the  candidate  solution. 

F  Superposition  can  be  used  for  adding  linear  differential 

equation  solutions. 

F  The  method  of  partial  fractions  is  used  to  solve  for  the 
particular  solution  of  a  differential  equation. 

F  The  number  "e"  is  a  variable. 

F  The  Laplace  transform  of  the  characteristic  equation  appears  i^ 
the  denominator  of  the  transfer  function.  m 

F  There  is  a  general  technique  vdiich  can  be  used  to  solve  any 
linear  differential  equation. 

F  Cramer's  rule  is  in  centimeters. 

F  Cramer's  rule  is  an  outdated  method  of  solving  simultaneous 

equations . 

F  The  treinsfer  function  conpletely  characterizes  a  linear  system. 

F  The  Heaviside  Expansion  theorem  is  often  cited  by  weight 

watchers. 

F  A  root  plot  is  a  short  hand  method  of  presenting  transient  time 
response . 

F  The  settling  time  is  a  measure  of  danping  ratio  of  a  system 

without  regard  for  the  danped  frequency. 

F  If  y  =  f(x),  then  y  is  the  dependent  and  x  the  independent 

vari^le . 
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3.32.  The  following  terms  are  inportant.  Define  and  provide  symbols  for 
those  you  are  not  sure  of. 

Differential  Equation 
Dependent  Variable 
Independent  Variable 
Ordinary  Differential  Equation 
Partial  Differential  Equation 
Exact  Differential  Equation 
Linear  Differential  Equation 
Degree  of  a  Differential  Equation 
Order  of  a  Differential  Equation 
General  solution 
Transient  solution 
Particular  solution 
Steady-state  solution 
Forcing  function 

Input  to  a  system  (related  to  the  differential  equation) 

Output  of  a  system  (related  to  the  differential  equation) 

Time  Constant 
Danping  ratio 
Danped  natural  frequency 
Natural  Frequency 
Undanped  response 
Underdamped  response 
Overdanped  response 
Unstable  system  response 
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Critical  Damping 

Linear  system 

Laplace  Transform 

Inverse  Laplace  Transform 

Itoit  Step 

IMit  Impulse 

Ramp  function 

Transfer  function 

Pole 

Zero 

Root  Plot 

Root  Locus 

Rise  time 

Settling  Time 

Peak  Overshoot 

Time  to  peak  overshoot 

Steady-state  error 

3.33.  Solve  the  following  problems.  Sketch  root  locus  plots,  and  find 
40^,  and  T  vdiere  appropriate. 


A. 


^y 

dx^ 


+  4x 


B.  dy 


-2xydx 
(x^  +  cos  y) 


C. 


dx 

at 


+  t^x 


0 
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D. 


E. 


F. 


G. 


dt^ 

dt^ 

^x 

dt^ 

dt^ 


-  5  ^  +  6x 
-4g.4x 

+  4x  =  0 

+  7  ^  +  22x 


0 


0 


0 


H.  Given:  =  2  sin  3x  +  2  cos  3x 

Find  A  and  <|>  in  the  expression 
=  A  sin(3x  +  <|>) 

The  following  problems  are  the  same  as  D  thru  G  with  forcing  functions. 


I. 


^x 

dt^ 


c  dx  .  , 

5  ^  +  6x 


=  9 


J. 


£x 

dt^ 


4^  .4X 


=  e 


-2t 


K. 


^X 

dt^ 


+  4x 


sin  3t 


L. 


^x 

dt^ 


+  7 


dx 

at 


+  22x 


t 
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The  following  problems  are  the  same  as  D  thru  G  with  forcing  function 
and  initial  conditions: 


M. 


N. 


0. 


dt^ 

6x  -  9 

o  o 

^x  _ 
dt^ 

■  ^  m 

+  4x  =  e^^ 

o  o 

+  4x  = 

sin  3t 

X(0) 

dt^ 

x(0) 

dt^ 

7  j. 

■  ^  at 

22x  =  t 

o  o 

3/2 

2 


2 

4 


3.34.  Solve  the  following  problems  using  Laplace  techniques: 


A.  - 

c  dx 

^  at 

+  6x  =  9 

x(0) 

=  3/2 

dt^ 

x(0) 

=  2 

B. 

+  4x  =  e^^ 

x(0) 

=  2 

dt^ 

x(0) 

=  4 

3.35.  Given  the  set 

of  equations 

t 


vhere  x(0)  =  y(0)  =  0,  find  y(t)  using  Laplace  transform  methods. 
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ANSWERS 


3.1, 

y 

3.2. 

y 

=  e-*  -  sin  «  X  ^  ^  , 

O) 

3.3. 

y 

X®  ''' 

=  336  -"  2  +C,*  +  C 

3.4. 

y^ 

=  -2x®  +  C 

3.5. 

=  Cx^  ^ 

3.6. 

Not 

exact. 

3.7. 

f  ^ 

x^ 

=  ^  +3xy  +  |  -y  +  C 

3.8. 

Not 

exact. 

3.9. 

f  = 

II 

+ 

to 

+ 

II 

O 

3.10, 

Yt 

=  Ce-‘^'* 

3.11. 

Yt 

^  ^  -3t  ^  8t 

=  C  +  C6  +C6 

3.12. 

Yt 

=  C^e"®^  +  Cjte”®’' 
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cos  (3t  + 


3.13.  Yt  =  C 

3.14.  y  =  -lle'^’^  +  12e"^^  -  3te"^’^ 

3 . 15 .  y  =  6  ~  •jg  to  +  cos  t  + 


3.16. 

y  = 

3e”*^  cos  3t  +  3/10  t 

3.17. 

y  = 

_e-2/3t 

+ 

3.18. 

VT 

C  = 

1.02 

3.19. 

w  = 

n 

2 

c  = 

1 

“d  = 

0 

3.20. 

“n  = 

VIo 

C  = 

0.316 

3.0 

3.21. 

X  = 

3/2 

6 

3.22. 

x(s) 

s^ 

+  36 

3s^ 

+  s  +  6 

3 

3.23. 

x(s) 

(s  +  1)^  +  9 

s^  -  2s  +  5 

sin  t 
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+  12s  +  1 


3.24. 


6  4s 

~2 - 2 - 2 

s^  (s^  +  4r 

4s^  +  3s  -  1 


3.25. 


2e‘^’=  -  3e"^^  +  6e'^ 


3.26. 


+  e“^  +  te"’' 


3.27.  y(t)  =  1  +  2/3  sin  3t  +  cos  N/s^  t  +  ^  sin  \/3/4  t 


3.28. 


=  26/5  e“  -  1/5  cos  t  +  2/5  sin  t 


3.29. 


-6e-^^  -  3te-^^  + 


3.30. 


1/4  (1  -  e~  (cos  2t  -  sin  2t)) 


=  1/4  (-1  e""'*  (cos  2t  +  3  sin  2t)) 


3.33.  A.  y  =  ^2  +  2/3  +  Cx  +  q 


B.  X  y  +  sin  y  =  C 


C.  X  =  Ce 


-tV4 


=  +  qe^‘ 
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F.  C  = 

0 

2 

«d  = 

2 

o 

• 

II 

4.69 

“d  = 

3.12 

0.746 

it 

< 

/8 

4*  = 

11/4 

I.  X  = 

q  +  Cj  +  3/2 

It 

x 

• 

+  1/2  t"e“ 

K.  X  = 

Cj  cos  2t,+  Cj  sin  2t  -  1/5  sin  3t 

■ 

L.  X  = 

(q  cos  t  +  Cj  sin-^  t)  + 

1/22  t  - 

M.  X  = 

-  2e“  +  2e^’'  +  3/2 

N.  X  = 

2e  +  1/2  t  e 

0.  X  = 

3/20  sin  2t  ^  1/5  sin  3t 

11 

X 

(-^  cos  1  +  .016  sin 

7 

484 

3.34.  A.  x(t)  =  I  - 

B.  x(t)  =  2e'“  +  1/2  t^e“ 

3.35.  y(t)  =  3t  - 

x(t)  =  1/2  -  t 
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